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In this paper, we develop a non-perturbation theory for describing decoherence dynamics of 
electron charges in a double quantum dot gated by electrodes. We extend the Feynman- Vernon 
influence functional theory to fermionic environments and derive an exact master equation for the 
reduced density matrix of electrons in the double dot for a general spectral density at arbitrary 
temperature and bias. We then investigate the decoherence dynamics of the double dot charge qubit 
with back-action of the reservoirs being fully taken into account. Time-dependent fluctuations and 
leakage effects induced from the dot-reservoir coupling are explicitly explored. The charge qubit 
dynamics from the Markovian to non-Markovian regime is systematically studied under various 
manipulating conditions. The decay behavior of charge qubit coherence and the corresponding 
relaxation time Ti and dephasing time T2 are analyzed in details. 
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I. INTRODUCTION 

Double quantum dot systems have been attracting 
much attention because of their intriguing properties and 
their potential applications in nanotcchnology and quan- 
tum information processing [J, [3]. The basic structure 

[ of a double quantum dot system can be viewed as elec- 
trons confined in an electrostatic potential of double wells 
created by the fabricated gates, source and drain elec- 
trodes in the heterostructure of a semiconductor. The 
heterostructure of GaAs/AlGaAs is a typical example 

] for the realization of gate-defined double dots. The tun- 
ability of various couplings and energy levels in the dots 

[ makes it a promising a quantum device (see, for exam- 
ples, [1, 0, H H, 0| ) • Maintaining electron coherence in a 
double quantum dot is an important ingredient in mak- 
ing it part of a quantum information processor. However, 
fluctuations and dissipations brought up by quantum op- 
erations of manipulations and measurements as well as 
various features of the material enrich the physics of dou- 

' ble dot systems more than perfect coherent evolution. 
Thus a lot of attentions have been paid to investigate 
how various noises and interactions with the surround- 
ings attenuate the coherent evolution of electrons in the 
double dot. 

In this paper, we will concentrate on a non- 
perturbative dynamical theory for charge qubit manip- 
ulation with a double dot system gated by electrodes. In 
the quantum computing scheme in terms of double dots 
where the electron charge degree of freedom is exploited, 
the effects in deviating the coherency of charge dynamics 
are summarized in the fluctuations of the inter-dot cou- 
pling and energy splitting between the two local charge 
states as well as the dissipation induced damping effects. 
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The amplitudes of these fluctuations can be estimated 
from measurements of the noise spectrum of electron cur- 
rents and the minimum line width of elastic current peak 
Parallel theoretical works have been developed with 
different approaches in the literature for the purposes 
of both simulating the experimental results and under- 
standing the physical mechanisms living in the double 
dot. In the present work, we shall extend the Feynman- 
Vernon influence functional theory Q to fermionic envi- 
ronments and derive an exact master equation describ- 
ing the coherent and decoherent dynamics of the electron 
charges in the double dot with the back-action effects of 
the reservoirs being fully taken into account. 

Stochastic noise processes resulted in a time depen- 
dent Hamiltonian for the double dot have been widely 
analyzed in simulating the charge dynamics under noise 
influences. The Bloch-type rate equations for describing 
the double dot transport properties have been investi- 
gated by Gurvitz and co-worker used the many-body 
Schrodinger equation approach which is further applied 
to study decoherence of double dot charge qubit in [l^l . 
The phonon assisted processes have been investigated 
within the Born-Markov regime by Brandes et at us- 
ing Born-Markov typed master equation [ll| . A general 
expression of the qubit density matrix in case of pure de- 
phasing was used by Fedichkin and Fedorov [l^l to 
study the error rate of the charge qubit. Stavrou and Hu 
[1^ considered in details the wavefunctions of the double 
dot charge qubit for decoherence analysis. Karrasch, et 
al. come with the functional renormalization group ap- 
proach in dealing with the transport aspects of the mul- 
tiple coupled dots (isj . The non-Markovian dynamics 
has also recently been studied by a suitable spin-boson 
model considering the acoustic phonons by Thorwart, et 
al. using numeric al q uasiadiabatic propagator path in- 
tegral scheme [l^, Il7l |. Without Born-Markov approxi- 
mation, Wu et al. had devised an analytical expression 
for the dynamical tunneling current using a perturbation 
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treatment based on a unitary transformation [l8|,ll9|. Ef- 
fects from Coulomb interaction between the dots and the 
gate electrodes with the formulation of kinetic equations 
have been presented by Woodford et al. [2^. The di- 
versity in the methodologies and issues concerned in the 
literature show the physical richness of this novel system. 

As we can sec there are many factors competing to 
play the consequent physics in the double quantum dot 
system. To single out one factor from the others on the 
resulted dynamical properties of this charge device, we 
shall first concentrate in this paper on the effects in- 
duced by dot-reservoir coupling, where the double dot 
is designed in the strong Coulomb blockade regime such 
that each dot only contains one energy level. The reser- 
voirs consist of the source and drain electrodes which are 
controllable through the bias voltage. A schematic plot 
of the system is shown in Fig. [T] Thus the total Hamil- 

I i I 

Source (jS^j) (Ootl) j^^g 
I ^ <3c G, ^ 

FIG. 1: A schematic plot for the double quantum dot system 
tonian of the system we concern in this work is given 

by 

H = Eia\ai + E2a\a2 + Tc{a\ai + a\a2) 
+ ^^{£Lka\^aLk + £Rka\^^aRk) 

k 

+ '^{tiLka{aLk + t2Rka'^f(k^2 + ii-c) (1) 

k 

which contains the Hamiltonians of the double dot, the 
source and drain electrodes plus the interaction (electron 
tunneling processes) between them. The notations follow 
the convention and will be specified in details later. 

Our treatment is based on the exact master equation 
we derived for a general spectral density of the electron 
reservoirs at arbitrary temperature and bias: 

+ (0(2ajpa| - ajajp - pajaj) 

where p is the reduced density matrix of the double dot 
obtained from the full density matrix of the double dot 
plus the reservoirs by tracing out the environmental de- 
grees of freedom, and 

H'{t) ^ E[{t)a\ai + E'2{t)ala2 + T'^{t)a\a2 + T';{t)a\ai 

(3) 



is the corresponding effective Hamiltonian. All the time- 
dependent coefficients in the above equations will be de- 
rived explicitly and non-perturbatively in Sec. IIF The 
time-dependent fluctuations of the energy levels, E'^{t), 
and the inter-dot transition amplitude, T^{t), are the 
rcnormalization effects risen from the electron tunnel- 
ing processes between the double dot and the reser- 
voirs. Other non-unitary terms describe the dissipa- 
tive and noise processes with time-dependent coefficients, 
F°(t),F(t) and TP{t), depicting the fuU non-Markovian 
decoherence dynamics. Eq. ^ is obtained without con- 
sidering the inter-dot Coulomb repulsion. But as we 
will show explicitly in Sec. Ill it is easy to extend to 
the strong inter-dot Coulomb repulsion regime where the 
strong inter-dot Coulomb repulsion simply leads one to 
exclude the states corresponding to a simultaneous occu- 
pation of two dots from Eq. ([2|). Then the exact master 
equation allows us to exploit the intrinsic quantum deco- 
herence effects in the electron charge coherency brought 
up by the tunneling processes between the dots and the 
reservoirs through the bias controls. 

The master equation ([2]) is derived by extending 
the Feynman- Vernon influence functional theory Q to 
fermion coherent state path integrals [2l[. Historically, 
since it was first developed by Feynman and Vernon in 
1963 for quantum Brownian motion (QBM) modelled as 
a central harmonic oscillator linearly coupled to a set of 
harmonic oscillators simulating the thermal bath, the in- 
fluence functional theory has been widely used to study 
dissipation dynamics in quantum tunneling problems |22l | 
and decoherence problems in quantum measurement the- 
ory [i^, [23|. In these early applications, the master equa- 
tion was derived for some particular class of ohmic envi- 
ronment [2^, [2^, [l^l . The exact master equation for the 
QBM with a general spectral density at arbitrary temper- 
ature was obtained by Hu and co-workers in 1992 [28| . 
Applications of the QBM exact master equation cover 
various topics, such as quantum decoherence, quantum- 
to-classical transition and quantum measurement theory, 
etc.ji^l- Very recently, such an exact master equation is 
further extended to the system of two entangled opti- 
cal fields and two entangled harmonic oscillators for the 
study of non-Markovian entanglement dynamics in quan- 
tum information processing [29l. [30|. Nevertheless, using 
the influence functional theory to derive the exact master 
equation has been largely focused on the bosonic type of 
environments, up to date. 

On the other hands, the development of quantum 
transport theory in nanosystems has continuously at- 
tracted attention in the last two decades because of the 
great achievements in nanotechnology, where the reser- 
voir is, in many cases, a fermion system. The traditional 
approach to study the quantum transport in nanosystems 
is Schwin ger- Keldysh's nonequilibrium Green function 
formulism [31I. [s^ which has been extensively used in suc- 
cessfully describing various quantum transport phenom- 
ena, such as Kondo effect, Fano resonance and Coulomb 
blockade effects in quantum dots [H, [H, [H, [s^ . Mas- 
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ter equations for quantum transport through quantum 
dots have also been derived, but mostly in t he p erturba- 
tion theory up to the second order [sj, [H, [3§, 113, [4l[ . 
The exact master equation can be, in principle, obtained 
through the real-time diagrammatic expansion approach 
developed by Schon and co-workers [H, S^], as shown 
recently by one of the authors and co-worker (43j . An- 
other interesting formulism is the recently published hi- 
erarchical expansion of the equations of motion for the 
reduced density matrix by Yan et. al. [44|. Nevertheless, 
in contrast to the bosonic environments |28| . an explicit 
formula of the exact master equation for fermionic en- 
vironments with a general spectral density at arbitrary 
temperature and bias has not been carried out except for 
Eq. ([2]) in this work. 

Unlike the quantum transport in nanosystems where 
people pay more attention on the tunneling current spec- 
trum and its statistics, one cares in quantum information 
processing how the qubit coherency can be maintained 
for fast quantum operations where a strong coupling is 
required. Then a nonperturbative (with respect to the 
coupling between the system and its environment) mas- 
ter equation is more desirable for the precision manip- 
ulations of qubit states. Eq. ([2]) obtained in this paper 
has fully taken into account the back-action of the elec- 
tron reservoirs at arbitrary temperature and bias. This 
master equation is also valid for a general spectral den- 
sity. It is fully non-perturbative that goes far beyond the 
Born-Markovian approximation often used in the litera- 
ture. It enables us to explore the dynamics of the electron 
charge coherence in the double dot, from Markovian to 
non-Markovian regime under various manipulating con- 
ditions. Many other approximated master equations that 
have been developed in the literature can be obtained at 
well defined limits of the present theory. 

The rest of the paper is organized as follows. In the 
next section, we will use fermion-coherent-state path in- 
tegral approach to solve exactly the electron dynamics in 
an isolated double quantum dots, as an illustration. We 
then extend the Feynman- Vernon influence functional 
theory originally built on the coordinate representation in 
quantum mechanics to fermion coherent state represen- 
tation. The exact master equation for the reduced den- 
sity matrix of the double dot system coupling to electron 
reservoirs is derived in Sec. Ill, where we also reproduce 
many other approximate master equations at well defined 
limits of the present formulae. In Sec. IV, we investigate 
the non-Markovian decoherence dynamics of this device 
including the tunneling induced fluctuations in the en- 
ergy splitting and inter-dot coupling of the double dot 
as well as the noise and dissipation effects on the charge 
mode populations and interferences. The leakage effect 
is also discussed together there. The decay behaviors of 
charge qubit coherence and the corresponding relaxation 
time Ti and dephasing time T2 are analyzed in details. 
Conclusive remarks are given in Sec. V, and Appendices 
are presented for some detailed derivations. 



II. FERMION COHERENT STATE PATH 
INTEGRAL APPROACH TO AN ISOLATED 
DOUBLE DOT 

To illustrate the fermion coherent state path integral 
approach to the electron dynamics in a double quantum 
dot, we shall consider in this section a simple solvable 
system, a single electron in an isolated double dot, before 
we go to explore the realistic system in the next section. 
We also assume that each of the dots contains only one 
energy level, Ei and E2 respectively. The Hamiltonian 
of this isolated double quantum dot is 

H = Eialai + E2a\a2 + Tc,{a\ai + a\a2). (4) 

The notations follow the convention: a\ 2 are the creation 
operators for electrons in the double dot, and Tc is the 
electron transition amplitude between the dots. In terms 
of the density operator, the time evolution of the system 
is described by 

p{t)^U{t-to)p{t^)U\t-t^), (5) 

where the density matrix p{t) is the state of the system 
at a later time t, p{to) is the state at the initial time to, 
U{t — to) = exp { — — io)} is the evolution operator 

of the system, and we let fi, = 1 hereafter. 

Using the fermion coherent state representation [i^ . 
the density matrix at time t is expressed as {C\p{t)\ri) = 
p{C*yVit) where ^ and 77 are two Grassmann variables 
characterizing the two-mode fermion coherent states, 

\0= n 1^^) ' I6>=exp(-e4)|0), (6) 

1=1,2 

The fermion coherent state defined above is an eigenstate 
of the fermion annihilation operator, a^j^i) = As 
these coherent states are over-complete, they obey the 
resolution of identity, / — I7 where the inte- 

gration measure is defined by dp{(,) = Y[i^~^^^^'^^i'^^i- 
Note that the fermionic coherent states we used here are 
not normalized, and the normalization factors are moved 
into the above integration measure. Moreover, these co- 
herent states are also nonorthogonal. The overlap of two 
fermionic coherent states is: = cxp(^t^') with a 

matrix notation = {^l ^2)- 

The use of the coherent-state representation makes the 
path integral formulation for a fermion system generally 
possible. In the fermionic coherent-state representation, 
the time evolution of the density matrix becomes 

P{Cf,Vf,t) = j dp(Co)dp{r]Q)p{^^,r]o,to) 

xJ(^;,77/,t|^,0,77*,to), (7) 

where JiCf^Vf^A^o^Vo^to) = {0\U{t - to)\Co){rio\UKt - 
to)\'i]f) is the propagating function in which the forward 
and backward transition amplitudes {S.f\U{t — to)\S,o) and 
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{r]o\W{t — to)\rif) can be solved exactly using the path 
integral for the Hamiltonian 

Explicitly, the fermion coherent state path integral for 
the forward transition amplitude is given by 

ioiuit-toMo) = / p[e,r]cxp(*5cr,e]), (s) 



where the action is 

i=l,2 



(9) 



In the above equation, the path integral in- 
tegrates over all pathes ^i(T) and £,*{t) bounded by 
^i{to) — ClO and £,*{t) ~ £*j with i ^ i' . Since the action 
in Eq. ([8]) has a quadratic form, the path integral can 
be exactly carried out with the stationary path method 
Slii. The result is 



exp 



El 



er/C«(i)+c(io)6o 



(10) 



where ^j*(io) and £,i{t) are determined by the solution of 
the equations of motion 

Ur) +i[E^Ur) + TS'ir)] = 0, (Ha) 
C(r) -i{E,^Ar) + T,i*\r)\ = 0, (lib) 

with the boundary conditions $i(io) = Cio and = 
where i ^ i' . The backward transition amplitude 
{i'j{)\U'^{t — t{))\rjf) can be found by the same procedure. 
From Eq. pT|) we can introduce a 2 x 2 matrix Uo(i) such 

that £,{t) = Uo{t)^o and ^''{to) = (,fUo{t) where Uo{t) 
satisfies the equation of motion 



uo{t) 



El T, 
Tc E2 



uo{t) = 



with the boundary condition {uo)ij{to) ~ 5ij. 
tion of IWi is: 



(cos (p{t)I — in - cr sin (^(t)) 



(12) 
The solu- 

(13) 



with 4){t) = E{t — to), ip{t) = flo{t — to)/2. Here we have 
also defined E = ■7){Ei + E2) and the Rabi frequency 
Hq = Ve^ + A^, where e = Ei — E2 is the energy level 
splitting in the double dot and A = 2Tc the inter-dot 
tunnel coupling, cr is the Pauli matrix, / a 2 x 2 identity 
matrix and the unit vector n = (A,0, e)/f2o. Then the 
propagating function becomes 



J{Cf,Vf,t\^o,Vo,to) = exp(^}uo(0^o + vWoii)^/) 



If we use further the D-algebra of fermion creation and 
annihilation operators in the fermion coherent state rep- 
resentation 1^) [2ll |. 

6ie)=a.|0, -^IO-«IlO, (15a) 



(15b) 



(14) 



it is easy to derive the equation of motion for the density 
operator p{t): 

ij 

One can see from Eq. IT^ that iio(t)wo(t)^^ = 
/ E T \ 

E^- )' "^^^^ simply reduces Eq. pB]) to the fa- 
miliar Liouvillian equation, 

p{t)^-^[H,p{t)], (17) 

as we expected. Having ensured the path integral tech- 
nique can reproduce the dynamical equation for an iso- 
lated double quantum dot, we will apply it to the double 
quantum dot coupling to electron reservoirs in the next 
section. 

III. THE MASTER EQUATION FOR A DOUBLE 
QUANTUM DOT GATED BY ELECTRODES 

A double quantum dot between two electron reservoirs, 
the source and drain electrodes controlled via a bias volt- 
age (see Fig.[T]), has a total Hamiltonian as 

H = H(iot + ^^lov + Hj, (18) 

where 

Hdot = Eia{ai + E2a\a2 4- Tc,{a\ai + a\a2) (19a) 
is the Hamiltonian of the double quantum dots, 

{£Lka\kO-Lk + £Rka^ii.k°'Rk) (19b) 

k 

is for the source and drain electrodes (reservoirs) and 

Hi = ^{tiLka\aLk + hRkaRk^-^ + H.c) (19c) 
fe 

is the coupling (interaction) between the double dot and 
the reservoirs that depicts the electron tunneling be- 
tween them, where subscript k labels an electron state 
in the reservoirs, L and R denote the reservoirs of the 
source (left) and drain (right) electrodes, respectively, 
and tiLj;{t2Rk) is the electron tunneling amplitude be- 
tween the source (drain) and the left (right) dot. Since we 
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only concern in this work the dynamics of charge qubit, 
we omitted spin degrees of freedom for electrons. 

It should be pointed out that in Eq. (|19ap we did not 
include explicitly the inter-dot Coulomb repulsion. This 
is because a typical inter-dot Coulomb energy is of the 
order of hundreds /^eV, which is much larger than the en- 
ergy level splitting e and the inter-dot tunnel coupling A 
(both are of the order of tens /ieV or less) for charge qubit 
manipulation [1]. As a result, the inter-dot Coulomb in- 
teraction simply leads one to exclude the states corre- 
sponding to a simultaneous occupation of the two dots 
S Thus, we will derive an exact master equation 
for the reduced density matrix of the double dot with- 
out considering the inter-dot Coulomb repulsion at be- 
ginning. The master equation of the double dot in the 
strong inter-dot Coulomb repulsion regime is then ob- 
tained from the exact master equation by explicitly ex- 
cluding the states of doubly occupied two dots in terms 
of Bloch-type rate equations, as we will see later. 

To derive non-puterbatively the master equation of the 
reduced density matrix for the double dot system, we 
adopt the treatment of the Feynman- Vernon influence 
functional theory Q. This approach within the frame- 
work of path integral traces over the degrees of freedom of 
the environment (here the electrodes) into a functional of 
the dynamical variables of the system (the double quan- 
tum dot). This functional is called the influence func- 
tional by Feynman and Vernon, which contains all the 
dynamical effects from the back-action of the reservoirs 
to the system due to the coupling between them. Fol- 
lowing the fermionic path integral technique presented 
in the last section, we depict the route to an exact mas- 
ter equation aided with the results derived in details in 
the appendices. 

A. Influence functional 

Explicitly, the total density matrix of the double dot 
plus the reservoirs obeys the quantum Liouvillian equa- 



tion idptot{t)/dt = [H, ptot{t)], which yields the formal 
solution: 

Ptot(i) = e-*^(*-*°Vtot(io)e^^(*-*''). (20) 

As we are interested in dynamics of the electrons in the 
double dot, we shall concentrate on the reduced density 
matrix for the double dot by tracing out the environ- 
mental variables: p{t) ~ tr£;ptot(t)- Assuming init ially 
the dots and the reservoirs are uncorrelated [13, |28| . 
Ptot(^o) = p{to) ^ pE(to), then the reduced density ma- 
trix describing the full dynamics of the electrons in the 
double dot becomes 

(OIpWI'?/) = J dp{Co)dpim){Mto)\m) 

^Jirf,Vf,t\^o,Vo,to), (21) 
where the propagating function is defined as 

J v[C^;v*vV^'^^^'' ^i-^*'"*^ ""^rC;^*^] (22) 

with Sc [C* 5 £,] being the action of the double dot in 
the fermion coherent state representation given by 
Eq. ([5]), and is the influence functional 

which takes fully into account the back-action ef- 
fects of the reservoirs to the double dot and modi- 
fies the original action of the system into an effec- 
tive one, e^('5=[«*^^l-^*[''*'''l)j^[rC;f/*?7] = e^-^="'K*«'''*''l. 
The path integral f [C*^;?7*^] integrates over all pathes 
£,* (t) , ^(t) , rj* {t) , and rjir) in the Grassmann space 
bounded by ^(0 = Cf, ^(to) = 6, V*{to) = Vo^ and 
?7(t) = r]f. 

Let the reservoirs be initially in thermal equilibrium 
states at temperature the influence functional can 
then be solved exactly with the result (see the derivation 
in the Appendix A): 



HC^;V*v]^cxpY^ {- fdr rdr'(i^a(r-r')C(r)e.(r') + i^l(r-r')?7*(T')77,(r)) 

+ fdr f dr' (Fa (r - r')< (r)6 (r') - (r - r') (r) + ^ (r)) (77, (r') + (r')) ) } • (23) 

J to J to 



The two time correlation functions in the influence func- tional, 

Fuir -r') = Y, |t,ifcpe---(^-^'), (24a) 

k 

Fi{T^T') = Y,f{ei,)\tuk\'e-''^'''^^-^\ (24b) 
fc 
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arc called the dissipation-fluctuation kernels, where 



with I = L, R for i = 1,2, respec- 



tively, are the fermi distribution functions of the electron 
reservoirs, and /ZL.i? are the corresponding chemical po- 
tentials. These nonlocal time-dependent functions con- 
tains the full dynamics effect from the reservoirs to the 
double dot. We should also point out that in the co- 
herent state representation, the influence functional for 
a many-fcrmion environment has a form similar to that 
of a many-boson environment [2^ except for some sign 
difference due to the antisymmetric properties of fcrmion 
degrees of freedom. 

The physical meaning of the above influence functional 
is very clear. The four terms contained in the exponential 
function of Eq. ((^51) correspond to four different physi- 
cal processes in a time-closed path formulism, see Fig. [2] 
The first term gives the contribution of the back-action 







■il3 

FIG. 2: The closed-time path for the trace of the density 
matrix of the reservoirs. 

effect of the reservoirs to the double dot in terms of the 
time correlation function Fii{T — t') of Eq. (|24ap in the 
forward process from the time to to the time t. For the 
double time integrals it contains, the first one starts from 
to to T and the second one from to to t, sums over all the 
time sequences of the these propagations. Resummation 
of these propagating processes up to all orders of the 
system-reservoir coupling results in exactly the exponen- 
tial function appearing in Eq. (|23[) . Similar to the first 
term, the second term are just the back-action effect of 
the reservoirs to the system in the backward process from 
the time t back to the time to which is just the complex 
conjugate of the first term. In terms of the Schwingcr- 
Keldysh Green-function approach, the time correlation 
function in these two propagating processes are just the 
time-ordered and anti-time ordered Green functions. The 
histories of the forward pathes ^* (r) , (r) and the back- 
ward pathes "Tj* [t) , rii{T) are mixed up through the time 
correlation functions as shown in the third and the fourth 
terms in Eq. ((23|) . The third term in Eq. (|23l) represents 



the mix of the forward and backward pathes at the time 
t. The last term is a mix of the forward and backward 
pathes at the time to where the initial equilibrium proper- 
ties of the reservoirs, i.e. the fermi statistics and the tem- 
perature of the reservoirs, naturally enter into the time 
correlation function F^i{t — t'), as shown by Eq. (|24b[) . 

As we see all the influences of the reservoirs on the 
double dot are embedded in the two time-correlation 
functions, i.e., the two dissipation-fluctuation kernels 
p4)) in the influence functional. These two dissipation- 
fluctuation kernels are related each other through the 
dissipation-fluctuation theorem. This will become clear 
by introducing a spectral density defined as 



E 

k 



\tuk\^5{uj - eik). 



(25) 



Obviously, the spectral density contains all the informa- 
tion about the reservoirs' density of states involved in the 
electron tunneling between the dots and reservoirs. Then 
the time-correlation functions (|24p can be expressed as 



J,;(^)e-'"("-"') 



(26a) 
(26b) 



This equation manifests the dissipation-fluctuation the- 
orem in an open quantum system. It tells that all the 
back-action effects of the reservoirs to the double dot are 
crucially determined by the spectral density Juiuj). 



B. Exact master equation 

Now we can derive the master equation for the re- 
duced density matrix. As wc sec the effective action 
after integrating out the environmental degrees of free- 
dom, i.e. combining Eqs. (22) and ((^ together, still has 
a quadratic form in terms of the dynamical variables of 
the fermion coherent states. Thus the path integral ([22]) 
can be solved exactly by utilizing the stationary path 
method and gaussian integrals [H, [i^. The resulting 
propagating function is simply given by 



J 



'/(e/,V/,i|Co,%*,to) = ^(i)exp ^ [ 



it) + e* (to)e.O , 77*o77, (to) + V* {tH f 



(27) 
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where A(t) is the contribution arisen from the fluctuations around the stationary pathes which will be given later. 
The stationary pathes ^i(i) and 77^(^0) ^-re determined by the equations of motion 



^,{r)+t[EMr)+m'{r)]+ dr' Fu{r - t%,{t') - / dr' Ff,{r - t')Ut') + ^,{t')] ^ 



r,,{T)+i[E,n,{T)+T,r,,,{T)]+ / dr' Fu{t - T')r^,{T') - / dr' Fu{t ^ t')[7j,{t') + ^,{t')] 



(28a) 



+ / dr'Fiir - r')[77,(r') + ^.(r')] = (28b) 

Jto 



subjected to the boundary condition ^i(to) = S,io and 
Vtit) = Vif, while er(^oj and Tj*{t) in Eq. ^ can 
be obtained from the complex conjugate equations of 
(j28bp and (|28ap . respectively, under the boundary con- 
dition ^*{t) = 4*/ ^iid ?7*(to) = ?7i'o- The local terms 
ii + iiE.^i + T^^v) and + i(i;,ry, + T^i^,,) in Eq. (US]) 
are intrinsic and well describe the coherent dynamics of 
the electron states in an isolated double quantum dot, 
as we have discussed in the last section. The nonlocal 
terms involving two different time correlation functions, 
Fii{T ~ t') and -fjf (t — t'), stem from the coupling to 
the reservoirs. These two time correlation functions play 
quite distinct roles in the equations of motion. The in- 
teraction between the double dot and the electron reser- 
voirs is mediated through electron tunnelings between 
them. The correlation Fii{T — t') describes the back- 
action of the reservoirs to the double dot due to the inter- 
action between them. However, the correlation function 
F^i{t — t') also harbors the Fermi-Dirac statistic effect 
of the electron reservoirs which exists even in the zero- 
temperature limit. The latter situation is quite different 
from a bosonic environment where F^iij — r') vanishes 
at zero temperature poj . 

The solutions of S,i{t), T]i{to) and 77* (t), £,*{to) deter- 
mined by Eq. ((28|) and its complex conjugate can be fac- 
torizcd from the corresponding boundary conditions. It 
is not too difflcult to find that 

m = (/ - v{t))-' [um„ + v{t)r,f] , (29a) 

v{to)^uHt)[I+{I-vit)r'vit)]rjf 

- [I - u\t){I - vit))-'uit)]^o- (29b) 

Here we have again used the matrix notation: = 
(Ci S.2) and a same form for 77. The new dynamical vari- 
ables expressed as two time-dependent 2x2 matrices, 
u(t) and v{t), satisfy the dissipation-fluctuation inte- 
grodifferential equations: 



(30a) 



u{t) + iMu{T) + / dr'G(r - r')u(r') = 0, 

Jto 

v{t) + iMv[T) + I dT'G{T -t')v{t') 



to 



dT'G'^{T~T')u{T'), (30b) 



with the boundary conditions Uij(<o) = Sij and (io) = 
0. while u{t) in Eq. (|30b|) obeys the backward equation 
of motion to u{t), namely, u{t) = u'' {t + to — r) for to < 
T < t. In Eq. we have also defined the 2x2 matrices 



E, 



) = Fu{t- 
- T')6ij. As we see w(t) 
r') while v{t) depends 
The full complexity of 



ij — j^iyc — i) and Tc{i ^ j 
T')%,and G^^{T^T')^Fi{T 
is determined purely by Fii{T - 
on both correlation functions, 
the non-Markovian dynamics of charge coherence in the 
double dot induced from its coupling to the reservoirs is 
thus manifested through these equations of motion. 

As £,*{to) and 77* (t) satisfying respectively the complex 
conjugate equations of 77; (t) and ^i(T) with the boundary 
condition £,*{t) = and ri*{to) ~ rjiQ, it is easy to find 
the similar solution from Eq. (^5]) for ^*{to) and r]*{t). 
Substituting these results into Eq. ((?7|) and note the fact 
that v{t) is a hermitian matrix at r = t, we obtain ex- 
plicitly the exact propagating function of the double dot 
system: 

J{Cf,Vf,t\^o, Vo,to) = A{t) exp {e} Ji(t)eo 



where 



Ji{t) = w{t)u{t), J2(i) = 
Jsit) = {t)w{t)u{t) - /, 



w{t) - I, (32a) 
A{t) = l/det(u;(i)). (32b) 



to 



with w{t) = (/ — v{t))^^. All these time-dependent coef- 
ficients can be fully determined by solving Eq. ([50)1 . 

Once the exact propagating function is obtained, the 
dynamics of the reduced density matrix, Eq. (j2ip . which 
takes fully into account the back-action of the electron 
reservoirs, can be completely solved for any given initial 
electron state of the double dot. The explicit solution of 
the reduced density matrix relies solely on the solution to 
the equations of motion (j30p which, in general, has to be 
solved numerically. To check the consistency, we may let 
the double dot be decoupled from the electron reservoirs, 
namely, set tuk = 0, then the dissipation-fluctuation ker- 
nels vanish. As a result, Eq. (|30aP is reduced to Eq. ([12]) 
and u{t) is reduced to uo{t) whose solution is given after 
Eq. ([T^ in the last section, while the solution of Eq. ()30bp 
gives v{t) = 0. Consequently the propagating function 
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Eq. ([31]) is reduced to Eq. ([T4|) which recovers the exact 
solution of the isolated double dot system shown in the 
last section. 

Having the explicit form of the propagating function 
pip in hands, it is straightforward to derive the mas- 
ter equation for the reduced density matrix directly from 
Eq. pT|) . Here we shaU deduce an operator form of the 
master equation such that aU the time-dependent coeffi- 
cients in the master equation are explicitly independent 
from the initial state of the double dot as well as from 
any specific representation. Taking the time derivative 
to Eq. (j3H) . chminating the initial state dependence and 
using the D-algebra of the fermion creation and annihila- 
tion operators in the fermion coherent state representa- 
tion, Eq. , the exact master equation of the double dot 
with a general spectral density at arbitrary temperature 
and bias is given by: 

p{t) =^{r!y(t)[ata^.,p(t)] 

+ Tij{t){2ajp{t)a\ - a\ajp{t) - p{t)a\aj) 

+ rfj- WKpWoI - a\p{t)aj - a\ajp(t) - p{t)a\aj)^ 

+ r"(t)pW, (33) 

where all time-dependent coefficients in Eqs. ([33]) are de- 
termined by u{t) and v{t) through the following relations: 

f^yW = ^["^"' - {u^Y^u%, (34a) 

r,,(i) = -\[uu-^ + {u^r^ii^j, (34b) 
rf^.(i) = [iiu-^v + w(ut)-iut _ y^^^^ (34c) 

and T^{t) = trF''. The first term in the master equation 
is indeed the generalized Liouvillian term which can be 
explicitly written as 

Y,n.,,{t)[ala,,p\^-i[H'{t),p] (35) 

where 

H'{t) = E[{t)a\ai + E'^{t)ala2 + T'^{t)a\a2 + (t)4ai 

(36) 

is an effective Hamiltonian of the double quantum dot 
with the shifted (renormalized) time-dependent energy 
levels and the shifted inter-dot transition amplitude, 

E[{t)^i^l,, (^-l,2), T'^{t)=inr2. (37a) 

Using the equation of motion (|30a[) . we further obtain 

E[{t)=E,-lxn[Wu{t)l 1 = 1,2, (37b) 

T'M = + '-[w^2{t) - w;^{t)i (37c) 



where 

(37d) 

It shows that the shifted energy levels, E[ 2(i), and the 
inter-dot transition amplitude, T'^[t), are entirely con- 
tributed by u{t) that involves only the time-correlation 
function Fu {t — t). The rest part in Eq. ([55|) describes 
the dissipation and noise processes of electron charges 
with non-Markovian behaviors having been embedded in 
these time-dependent 2x2 matrix coefficients, T(t) and 
V-^it). We caU T[t) and r^(t) the dissipation-fluctuation 
matrix coefficients or simply the dissipation-fluctuation 
coefficients hereafter. Note that T{t) is solely determined 
by u(i), and TP{t) is given by both u{t) and v{t). Thus 
all the time-dependent coefficients in the master equa- 
tion are non-perturbativcly determined by Eq. (|30p and 
fully account the back-action effects of the reservoirs to 
the double dot. 



C. Bloch-type Rate equations for zero and strong 
inter-dot Coulomb repulsion cases and the 
corresponding Markovian limits 

(a) No inter-dot Coulomb repulsion: To closely exam- 
ine the decoherence of electron charge dynamics in the 
double dot, it is more convenient to rewrite the master 
equation ([33)) as a set of Bloch-type rate equations in 
terms of the localized charge states in the double dot. 
Without considering the inter-dot Coulomb repulsion, 
the rate equations can be obtained directly from the 
master equation ([55)) in the charge configuration space 
containing the states of empty double dot, the first dot 
occupied, the second dot occupied and both dots occu- 
pied. We label these four states by = 0, 1,2,3, re- 
spectively. Then the master equation in the above basis 
becomes 

Poo = TllPll + r2ipi2 + ri2/921 + r22P22 + T^POO, 

(38a) 

pii = (r° - 2f ii)pii -I- sipi2 + s_p2i 

- rf'iPoO + f22P33, (38b) 
P22 = (r° - 2f 22)P22 + '^*+Pl2 + S+P2I 

- r22P00 + f 11P33, (38c) 
P12 = [-«£' + r° - trf]pi2 + S+pii + S_p22 

- TiaPoo + f 12P33, (38d) 

P33 = 1 12P2I + 1 21P12 - 1 I1P22 - 1 22P1I 

+ (po - 2trf )p33. (38e) 

Here the density matrix elements are defined by pij ~ 
{i\p\i)- We have also defined all the time-dependent co- 
efficients in the rate equations as V{t) = 2T{t) + r'^(i). 
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f (i) = r{t) + r^{t), 2± = ±iT,'(t) - f 12 and e'{t) = 
E[{t) - E'^{t) where E[(t), T^{t), T{t),Tf^{t) and T°{t) arc 
the time-dependent transport coefRcients contained in 
the master equation and are expHcitly given by Eqs. (j34p 
and ((37)) . These time-dependent coefBcients in the rate 
equations thus fuUy characterize the non-Markovian dy- 
namics of electrons in the double dot. The first and the 
last rate equations account electron charge leakage effects 
in this device, while other three rate equations depict 
charge qubit decoherence dynamics under the influence 
of the reservoirs. 

For a constant spectral density that has been widely 
used in the literature, the time-correlation function be- 



comes 



' 271 



doje 



o{t-t') 



r,5(T-r'), 



(39) 



where T; = 2TTpi\tak\^, {I = L,R for i = 1,2), and pl.r 
arc the densities of states for the left and right electron 
reservoirs. Then in the Markovian limit [t — t^, t — t^ ^ 
oo) [131, the integral kernels in Eq. ([50)1 reduce to 

dT'Fa{T - r')My(r') ^ ]-Tiu,j{t), (40a) 

tn ^ 

dT'^;^(r - T')u,,{r') ^ fiViu,,{T). (40b) 

to 

The factor ^ in Eq. (|40ap comes from the boundary of the 
time integration sitting upon r, and /l._r in Eq. (|40bp 
are the fermi distribution functions of the reservoirs. In 
this Markovian limit, the equation of motion (|30|) can be 
solved analytically and the solution is 



u(t) = exp < — 



iTr 



lE2 + ^ 



v{t) 



L 



(41a) 

u(r' + io)- 
(41b) 



With the above explicit solution for u{t) and v{t), we 
can obtain analytically all the coefficients in the master 
equation as well. Note that u{t) = u^t + to — t) depend- 
ing explicitly on the time t implies that dTv{T)\T=t 7^ 
dtv{t) = v{t), it is easy to find 



iiu ^ = — 



iE 



1+2 



iiu + v{v)) ^v) 



2+ 2 







IrTr 



(42a) 
(42b) 



As a result, for a constant spectral density in the Marko- 
vian limit, all the time-dependent coefHcients in the mas- 
ter equation become constant: 



(43a) 
(43b) 



In other words, in the Markovian limit with a constant 
spectral density, there is no renormalization effect to 
the energy level shift and the inter-dot transition ampli- 
tude. The dissipation-fluctuation effects are simply re- 
duced to the time-independent tunneling rates between 
the dots and the reservoirs. In fact, the differences of 
time-dependent and time-independent coefficients mani- 
fest the non-Markovian dynamics of electron charges in 
the double dot. We shall present quantitatively such dif- 
ferences in the next section. 

For a constant spectral density in the Markovian limit, 
the rate equations are simply reduced to 

Poo Il^lPii + Ib^rP22 - {fh^L + IrJ^r^Poo, (44a) 
Pii = -{Jl^l + /_Rr_R)pii 4- iTcXpu - P21) 

+ /LriP00 + /flrfl/933, (44b) 
P22 = -{IrXr + fL^L)p22 - iTc{pi2 - P21) 

+ fR^RPoo + IlTlPss , (44c) 
Pi2 = (-«£ ^ -)pi2 + iTcipii - P22), (44d) 

P33 = -{IiXl + /flFflJ/Oas -t- fL^LP22 + Ir^RPiI, 

(44e) 

where Jl.r = (1 ^ /l,/?) and e = Ei — E2- Under the 
large bias limit, Jl ==1, /r — 0, the above rate equa- 
tions reproduce the rate equations obtained by Gurvitz 
and Prager Q for the double dot without considering the 
inter-dot Coulomb repulsion: 



Poo = — FlPoo + ^RP22, 

Pii = TlPoo + ^rP33 + iTc{pi2 - P21), 

P22 = -(Fl + ri^)p22 - iTc{pi2 - P21), 



P12 = {-is 



2 



-)pi2 + iTc{pii - P22), 



P33 



~rflP33 + FlP22- 



(45a) 
(45b) 
(45c) 

(45d) 
(45e) 



(b) Strong inter- dot Coulomb repulsion: On the other 
hand, realistic experiments of the double dot are set up 
in the strong Coulomb blockade regime where not only 
each of dots has only one effective energy level but also 
there is no states of simultaneous occupation of the two 
dots. In other words, the configuration space of the local- 
ized charge states in the double dot system with a strong 
inter-dot Coulomb repulsion only contains the states of 
empty double dot, the first dot occupied and the sec- 
ond dot occupied, denoted by = 0,1,2, respec- 
tively. The corresponding rate equations in this strong 
Coulomb blockade regime for an arbitrary spectral den- 
sity can also be obtained by simply excluding the doubly 
occupied states in Eq. (|38p . Since the rate of the dou- 
bly occupied state |3) in the case of ignoring inter-dot 
Coulomb repulsion depends on the populations pn and 
P22- This probability flow from the states |1) and |2) to 
|3) should be redirected back into the states |1) and |2) 
in the strong inter-dot Coulomb repulsion regime. Mean- 
while, to ensure the probability conservation without the 
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doubly occupied states, a correction to the dependence of 
the coherence elements pi2 and p2i hi the rate equations 
for pii and p22 raust also be taken into account guided 
by the condition: 

rf2P2i + r^iPi2 - rfiP22 - r^2Pii = O- (46) 

This condition indeed forces the doubly occupied state 
to decouple from other states in the double dot, as one 
can see from the rate equation P33 in Eq. ([55)1 . These 
modifications can be done explicitly by taking the fol- 
lowing shift to the coefficients in the non-interacting rate 
equations fu Tu + ^rj, ^ - ipf^. with 

i j. In fact, the above coefficient shift also automati- 
cally cancels the r° dependence of P33 in Eq. , which 
is indeed a criterion for entirely excluding the double oc- 
cupied state from the reduced density matrix, as one can 
directly see from the expression of the master equation 
([55)1 . Then the rate equation of in Eq. ([55)1 becomes 
P33 = (— 2trr)p33, its solution is pasit) = if initially 
P33{to) = 0. This impHes that no leakage into the double 
occupied state will occur. As a result, the rate equa- 
tions in the strong inter-dot Coulomb repulsion regime 
arc given by 

Poo = Tiipii -f- T21P12 + ri2P21 + r22P22 + T'^POO, 

(47a) 

Pii = -f iiPii + Slpi2 + S_/92i - rfiPoo, (47b) 

P22 = -f 22/022 + '^+Pl2 + S+P2I - r22P00, (47c) 

P12 = [-ie' - ^trf]pi2 + S+pii + S_p22 - ^uPoo, 

(47d) 

where S±(t) = ±iT^{t) - i2(<). This set of the rate 
equations depicts the full non-Markovian dynamics of 
the double dot in the strong inter-dot Coulomb repul- 
sion regime. 

For the case of a constant spectral density in the 
Markovian limit, fjj(i) = {2r{t) +T^{t))ij fiTiS^j, 
S±(t) — *■ ±iTc and e'{t) e. The rate equations under 
such circumstances are reduced to the rate equations for 
the double dot in the strong inter-dot Coulomb repulsion 
regime, given in [ll[, 

Poo = -(/lFl + /j?r_R)poO + JiXlPw + fR^RP22, 

(48a) 

Pii = /lFlPoo + iTc{pi2 - P21) - /lTlph, (48b) 

P22 = fRTRPm ~ iTc{pi2 - P21) - fR^RP22, (48c) 

P12 = [-ie- ^(/lFl + ./i?rfl)]pi2 + iTcipii - ^22)- 

(48d) 

Furthermore, in the large bias limit, = 1, = 0, the 
above rate equations lead to Stoof-Nazarov's rate equa- 



tions Hi]: 

Poo = -FlPoO + rij;p22, (49a) 

Pii = FlPoo + iTcXpi2 - P21), (49b) 
P22 = -rfl,p22 - iTc{pi2 - P21), (49c) 

p 

P12 = (-«£ ^)pi2 + iTcipii - P22)- (49d) 

In a summary, we have derived in this section an ex- 
act master equation for the double dot gated by elec- 
trodes, and the corresponding Bloch-type rate equations 
psp without considering the inter-dot Coulomb interac- 
tion as well as the rate equations (|47)) for the strong 
intcr-dot Coulomb repulsion. For convenience, we call 
Eq. ([55)) the interaction-free rate equations and Eq. ([T7|) 
the strong-interaction rate equations hereafter. Other 
approximated rate equations that have been used in the 
literature are obtained at well defined limit of the present 
formulae. 



IV. NON-MARKOVIAN DYNAMICS OF 
CHARGE QUBIT 

With the above formulae we can now systematically 
explore the non-Markovian dynamics of the charge qubit 
for this double dot system. The coherence (decoherence) 
dynamics of electron charges in the double dot is deter- 
mined by its internal structure as well as external opera- 
tions. The internal structure includes the spectral prop- 
erties of the reservoirs as well as the couplings between 
the dots and the reservoirs embedded in the spectral den- 
sity. The external operations include charge qubit initial- 
ization, its coherence manipulation and the qubit state 
readout through the bias controls of the source and drain 
electrodes. Non-Markovian decoherence effects of these 
internal structure and external operations to the charge 
qubit are manifested through the time-dependent coef- 
ficients in the master equation, which is completely de- 
termined by Eq. ([50)1 after the spectral density Jii{u)) is 
specified. In fact, the time correlation functions directly 
tell us the length of the correlation time which determines 
to what extent the time-dependent fluctuation and mem- 
ory effect become important. The longer the correlation 
time is, the more memory effect acts on the electron dy- 
namics in the double dot and vice versa. 

To be more specific, we should first specify the spectral 
density Ja{uj) for the source and drain electrodes. Unlike 
the bosonic environment where a general spectral density 

J{lo) = riuj\^\ e (tjj, is a high frequency cutoff 

and 77 is a dimensionless coupling constant) was defined 
and used to classify the bosonic environment as Ohmic 
if n = 1, sub-Ohmic if < n < 1, and super-Ohmic 
for n > 1 [2^ . for a fermionic environment a general 
spectral density should not be a Poisson or Gaussian type 
distribution function because of the fermi statistics. Here 
we shall use a Lorentzian spectral density that has been 
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used in the study of the influence of a measuring lead on 
a single dot [491 and molecular wires coupling to electron 
reservoirs [50|. The Lorentzian spectral density wc used 
here has a form: 



Tidf /2Tr 



(50) 



where Ei is chosen to be the energy levels of the dou- 
ble dot. and / = L, i? for i = 1, 2, respectively. There 
arc two parameters in Jii{to) that characterize the time 
scales of the reservoirs. The parameters cZl,/? describe 
the widths of the Lorentzian distributions, which tell how 
many states in the reservoirs around Ei^2 effectively in- 
volve in the electron tunneling between the reservoirs and 
dots. Hence, the inverse characterize the time scales 
of the source and drain electrodes. Another parameter 
is the electron tunneling strength or the tunneling rate 
between the reservoirs and dots, T]^ fj, its inverse char- 
acterizes the time scale of the electron tunneling process 
itself between the reservoirs and dots. Indeed, Tl^h also 
describe leakage effects of electrons from dots to the reser- 
voirs and vice versa. Thus a Lorentzian spectral density 
well depict the time scales of non-Markovian processes in 
this open quantum system. 

Also, the choice of a Lorentzian spectral density makes 
it easy to recover the constant spectral density which has 
been often used in the literature. In fact, taking the large 
width limit, namely assuming all the electron states in 
the reservoirs has an equal possibility for electron tunnel- 
ing between the reservoir and dot, then 2TrJii{Lu) F; 
reproduces the constant spectral density that has often 
been used in the study of both quantum transport and 
quantum decoherence phenomena in nanostructures. In 
this limit, the time scale of the reservoirs is suppressed. 
Furthermore, in the previous investigations, especially 
in the study of quantum transport phenomena, one also 
takes a long time limit. Combining these two limits (the 
constant spectral density and long time limit) together, 
the exact master equation is reduced to the Bloch-type 
rate equations in the Markovian approximation obtained 
by others [1, [HI, HI] , as we have shown in the last section. 
Hence, with a Lorentzian spectral density, it is not only 
convenient to analyze in details the non-Markovian dy- 
namics but also enables us to easily make a comparison 
with the Markovian dynamics. 

For a Lorentzian spectral density, the corresponding 
temperature-independent time correlation function can 
be exactly calculated. The result is 



Fu{t 



^^exp{-(rf,+zi?,)|r-r'|}. (51) 



Obviously, rf^ describes the correlation times of the 
reservoirs. The wider/narrower it is, the shorter/longer 
the correlation time will be. The internal structure of the 
double dot is characterized by the energy level splitting 
e ~ El — E2 and the inter-dot tunnel coupling A = 2Tc, 
its time scale Tq is given by the inverse of the bare 



Rabi frequency f^o = V s'^ + A^. The non-Markovian dy- 
namics should be dominated when the two typical time 
scales, ^ ^^'^ '^l i? ' ^^'^ ™ ^^i^ same order of magnitude. 
There is another time scale, the reservoirs' temperature 
(3=1/ kT that also influences the non-Markovian dynam- 
ics of the charge qubit in certain cases. In the current 
experiments for charge qubit manipulation [J, the tem- 
perature is roughly fixed at 100 mK. We will take this 
temperature throughout our analysis to the charge qubit 
decoherence. 



A. Time dependent coefficients in the master 
equation and non-Marl<ovian dynamics 

Once the spectral density is specified, the full non- 
Markovian dynamics of charge qubit in the double dot 
can be depicted using the master equation (|33p. or more 
specifically the corresponding Bloch-type rate equations 
([38)) and ([47[) for the cases of no inter-dot Coulomb re- 
pulsion and strong inter-dot Coulomb repulsion double 
dot, respectively. To solve the master equation or equiv- 
alently the rate equations, we must determine first the 
time-dependent coefficients contained in these equations, 
namely the shifted (renormalized) energy level splitting 
e'{t) = E[{t) — E2{t) and inter-dot tunneling couphng 
A'(t) = 2r^(t), as well as the dissipation-fluctuation co- 
efficients F(t) and r'^(i). These transport coefficients are 
completely determined by the functions uit) and v{t) as 
the solutions of the dissipation-fluctuation equations of 
motion ([5(7[l which have to be solved numerically for a 
given spectral density. 

Using the Lorentzian spectral density ([50[) . wc can cal- 
culate explicitly all the time-dependent transport coeffi- 
cients in the master equations and then discuss the cor- 
responding non-Markovian dynamics by comparing with 
the Markovian limit in various different time scales. We 
shall first analyze the time-dependent coefficients for the 
charge qubit initialization where a bias is applied to the 
double dot and the double dot is adjusted to be off- 
resonance, i.e. eVsD = Ml — P^r 7^ and e = Ei — E2 ^ 
Q . We will examine when the large bias limit is reached 
and how the initialization works. After that we will go to 
the coherent manipulation regime where the double dot is 
set up symmetrically {Ei = E2) and the chemical poten- 
tials of the electron reservoirs are aligned above the en- 
ergy levels of two dots with zero bias voltage (/i^ = Mi?)- 
The time dependence of e'{t),T^{t) as well as T{t) and 
F'^(t) in this regime will tell us when the non-Markovian 
dynamics becomes important during the charge qubit 
evolution. 

The dissipation-fluctuation equations of motion (j30p 
show that only the solution of v{t) depends on the fermi 
distribution function in the reservoirs. In other words, 
only r^{t) sensitively depends on the bias. Other co- 
efficients, e'{t),T^{t) and T{t), do not depend on the 
chemical potentials hl,r and thus the bias. Their time 
dependencies are completely determined by the inter- 
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nal parameters of the double dot and the spectral den- 
sity. In the initialization scheme where a bias is pre- 
sented and the energy splitting of the two levels in the 
double dot is nonzero within the transport window, the 
time-dependence of is plotted in Fig. [3] by varying 
the bias voltage. We find that the large bias limit is 
reached at about 100 ficY for the given internal parame- 
ters e = 30 AicV, A = 10/^eV, Tl^r = A and dL,R = A/2, 
with JIq = Ve^ + A^ ~ 32/^eV being the bare Rabi fre- 
quency of the charge qubit. This large bias limit is shown 
in Fig. [3] that the curves for eVso = 90 fj.eV are very 
close to the curves of eVsD = 430 ^eV, and the curves 
for eVsD = 130 iJ,eV perfectly overlap with the curves of 
eVsD = 430 iieV. 




eVs 




— — - eVs 


D=90neV 


eVg 


D=130^eV 


— - - ■ eVg 


D=430).ieV 




FIG. 3: The time-dependence of V'^ by varying the bias 
voltage. The double dot is set to be off-resonance with 
e = 30 AteV, A = lO/xeV, Fl.r = A and dL,R = Oo/2. (a) and 
(d) show the two diagonal matrix elements of P''. Since F*^ 
is hermitian only one of the off-diagonal elements is plotted 
with the real part shown in (b) and the imaginary part in (c) . 
The period of the bare Rabi cycle is To ~ 2n/Qo. 



Fig. 3] shows the time-dependence of other coefficients 
for different energy level splitting of the double dot 
states (e = 0, 10, 70/ieV) with different tunneling rate 
(Fi = Tr = 2.5, 10,25/^eV). The result shows that 
e'{t) = for e = 0, and the initial value Re(T^')(0) al- 
ways equals to Tc with lm(r^')(0) being zero in all the 
cases we have calculated. For a small tunneling rate 
(^L,R < A/2), the time-dependence of these coefficients 
is almost negligible (see red solid lines in Fig. 5]). The 
time-dependent effect appears when the tunneling rates 
between the reservoirs and dots become relatively large. 
These time-dependent behaviors vary sensitively on the 
energy splitting of the double dot states. Increasing the 
energy splitting e changes the time-dependent behaviors 
of all the coefficients significantly, as shown in Fig 2] 
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FIG. 4: The time-dependence of the coefficients e'{t),T'^{t) 
and T(t) for different energy splitting e and tunneling rate 
Fl,r at the inter-dot coupling A = lOpeV and the spectral 
widthes dL,R. ~ A/2. 



The dynamics of electron charges in this initialization 
regime is plotted in Fig. [51 If we dope only one excess 
electron in the left dot with the right dot being empty, 
the efficiency of keeping this initial state decreases as the 
bias decreasing [see Fig. [5{a)]. If the two levels of the 
double dot are in resonance (e = 0) or the inter-dot tun- 
nel coupling is larger than the tunneling rates A > 2Tl^r, 
it also has a low efficiency to keep the double dot in the 
initial state pu = 1. A large bias configuration {eVsD 
larger than 100 /ieV) maintains the double dot in the 
initial state pn = 1 very well. Furthermore, for a large 
bias it also quickly leads the double dot into the state 
Pii ~ 1 even if the initial state is poo = 1 (both dots 
arc empty initially) or p22 = l(thc left dot is empty 
but the right dot is occupied by one excess electron), 
as shown in Fig. [5Kb). These numerical solutions are 
obtained using the interaction- free rate equations (|38p . 
But the strong-interaction rate equations (|47[) give qual- 
itatively the same result for initialization. Taking the 
bias to be 650/ieV and the reservoirs' ternperature to be 
lOOniK that have used in experiments [^, a very effi- 
cient initialization of the charge qubit can be obtained, 
as shown in Fig. [5] 

Meanwhile, Fig.[3|shows that for a relatively large tun- 
neling rate Tl^r, the smooth time oscillation of all the 
coefficients at small e become discontinuing at a relatively 
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FIG. 5: (a) The time evolution of pn, the probability of find- 
ing one excess electron in the left dot with the right dot be- 
ing empty, is plotted at various bias when the initial state is 
pii = 1. (b) The time evolutions of pn (black dash-dotted 
line) and poo (green long dashed line) when the initial state 
is an empty state (poo = 1), where eVsD = 650peV. Other 
input parameters used here are the same as that in Fig. [3] 



large e value. Such discontinuities correspond to the elec- 
tron hopping to the localized charge states in the double 
dot where the inter-dot tunnel coupling A is almost neg- 
ligible in comparison with the level splitting e. Such dis- 
continuities are also manifested perfectly in the electron 
dynamics. The discontinuities coincide with the times at 
which the electron is found in a localized charge state 
of the double dot in a very high probability. Fig. [DJa)- 
(b) are obtained using the rate equations ([55)1 and P7)) . 
respectively, where we plot the corresponding electron 
charge dynamics together with the time-dependence of 
Fii. We find that for a large bias, the electron charge dy- 
namics given by the interaction-free rate equations ([35]) 
and the strong-interaction rate equations (|47)) display the 
same feature, including the coincidence between the dis- 
continuities in the tinic-dcpcndcnt coefficients and the 
emergence of a localized charge state in the double dot. 
Note that if a sufficiently large bias is applied across the 
double dot, initialization can still be achieved regardless 
of these discontinuities. This is because a large bias is the 
most dominant factor in this situation. As a conclusion, 
initialization of charge qubit in the double dot can be 
easily achieved in a large bias limit with relatively small 
tunnel coupling and tunneling rates. All these parame- 
ters are tunable in experiments. 

Now we turn into the regime for charge qubit rota- 
tions where the double dot is set up at the resonant 
levels (El ^ E2 ~ E) and the fermi surfaces of the 
electron reservoirs are aligned above the resonant lev- 
els = fJ-R = M s-nd fi — E > 0). In other words, 
the double dot is set to be symmetric and unbiased for 
charge coherence manipulation Q. Figs. [7][n] show the 
time-dependencies of the shifted inter-dot tunnel cou- 
pling A'{t), as well as the dissipation- fluctuation matrices 
r(i) and T^{t) for the symmetric double dot by vary- 
ing the chemical potential fi — E, the spectral widths 
dL.R and the tunneling rates T^^n, respectively, from 
which we can determine the time scales within which 
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FIG. 6: The coincidence between the discontinuities in rii(t) 
(the green dot-dashed line) and the emergence of a localized 
charge state in the double dot. The input parameters e = 
70^eV, A = lOpeV, dL,R = A/2 and Pi, a = 25peV with 
the bias eVsD = 120^eV such that the initialization can be 
achieved, (a) for the double dot without considering the inter- 
dot Coulomb repulsion, (b) for the double dot in the strong 
inter-dot Coulomb repulsion regime. 



non-Markovian processes dominate the charge coherence 
dynamics at zero bias. We find that for the symmetric 
double dot, the shifted energy level splitting e'{t) and the 
imaginary part of the shifted inter-dot tunnel coupling 
ImA'(i) keep to be zero as we have already pointed out 
in Fig. m The off-diagonal element Fi2(t) and the imag- 
inary part of F^2(0 

are found also to be zero, while the 



diagonal elements Tii{t) = F22(i) and ^^^^{t) — F22(i)- 
All the time-dependent coefficients change in time in the 
beginning and then approach to an asymptotic value (the 
Markovian limit) at different time scales. These time- 
dependence behaviors will be used to analyze the deco- 
herence dynamics of charge qubit in the next subsection. 

In Fig. [7] we plot the shifted inter-dot tunnel cou- 
pling A'{t), the dissipation-fluctuation matrices Tii{t) 
and F^j^ (0 by varying the chemical potentials with re- 
spect to the energy levels of the double dot, fi~ E. The 
red solid, blue long dashed and the black short dashed 
lines correspond to /i— E' = 0, 25 and 50/ieV respectively. 
In Fig.[7]the spectral widths dL,R = where rig = A 

is the bare Rabi frequency of the double dot. In other 
words, the time scale of the reservoirs is chosen about two 
Rabi cycles of the system. A large amplitude variation 
of these time-dependent coefficients within the time scale 
of the reservoirs is clearly shown up in the figure. After 
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that time, all these time-dependent coefhcients approach 
to a steady value which corresponds to their asymptotic 
values as a Markovian limit. This indicates that the pos- 
sible non-Markovian dynamics is mainly caused by the 
time-fluctuations of these coefficients within the charac- 
teristic time of the reservoirs. 
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FIG. 7: (a) the shifted inter-dot tunnel coupling A'{t), (b)- 
(d) the dissipation-fluctuation matrix elements Tii{t), and 
rfi i2(^) by varying the chemical potential jj, — ij.l,r- The 
dot parameters E = Ei = E2 and A = 10 /ieV. Tl,r = A 
and dL,R = A/2. The matrix is hermitian so that only one 
of the off-diagonal elements is presented here. The red solid, 
blue long dashed and the black short dashed lines correspond 
to /i — iJ = 0, 25 and 50 neV respectively. 



Fig. [5] is the same plot with varying the electron tun- 
neling rates between the reservoirs and dots but fixing 
the chemical potentials at yU — i? = A. The small tun- 
neling rate Tl^h = 2.5 fj.eV (< A/2) does not show a 
significant time variation in the dissipation-fiuctuation 
coefficients, as also shown in Fig. 01 The shifted inter- 
dot tunnel coupling at this tunneling rate is very close to 
the bare one ( ~ A), sec the red solid line in Fig. [Hta). 
This implies that a small electron tunneling rate between 
the reservoirs and dots (corresponds to a small leak- 
age effect) does not manifest the non-Markovian dynam- 
ics significantly. Increasing the tunneling rates enlarges 
the charge leakage effect from the reservoirs to dots and 
vice versa, thus enhances the non-Markovian effects as 
well, as shown by the giggling and wiggling time evo- 
lutions of the dissipation-fiuctuation coefficients in the 
figure. The time-dependence of shifted intcr-dot tunnel 
coupling A' at large tunneling rates also has a signifi- 
cant shift from the bare one A besides the oscillation 
within the non-Markovian time region. Meanwhile, the 
charge oscillation frequency has different shifts from the 
bare Rabi frequency for different Tl^h values. How these 
time-dependent (non-Markovian) behaviors influent the 
charge coherence will be discussed in details in the next 
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FIG. 8: The same plot as in Fig. [7] with flxing the chemical 
potentials at /x — -B = A and varying the electron tunnel- 
ing rates between the reservoirs and dots Tl,r symmetrically. 
The other parameters are the same as in Fig. [T] 



Fig. [S] shows how the time-dependence of the shifted 
inter-dot tunnel coupling and the dissipation-fiuctuation 
coefficients change by varying the spectral widths cIl.r- 
We plot these time-dependent coefficients for three dif- 
ferent spectral widths: dL,R = 1, 5 and 25 /icV. The 
result shows that when dL,R = 25 /icV, the dynamics of 
the double dot already reaches to the Markovian limit, 
namely all the time-dependent coefficients approach to 
their asymptotic values in a very short time (less than a 
half cycle of the bare Rabi oscillation). When the spectral 
widths (Il.r. becomes small (< A) so that the character- 
istic time of the reservoirs becomes long), the time os- 
cillation of all the transport coefficients becomes strong. 
Correspondingly the charge dynamics is dominated by 
non-Markovian processes. 

The above numerical results tell that the spectral 
widths d^^R of the reservoirs (mainly as a memory effect) 
and the tunneling rates T^^r between the reservoirs and 
dots (mainly as a leakage effect) are two basic parameters 
to characterize the occurrence of non-Markovian dynam- 
ics in this double dot device. Comparing the results of 
Figs. [71 [51 and [HI we find that for the coherence manip- 
ulation of charge qubit where the the double dot is un- 
biased [3, the time for the dissipation-fiuctuation coeffi- 
cients reaching to a steady limit depends on the spectral 
widths of the tunneling spectra as well as the tunneling 
rates between the reservoirs and dots. The Markovian 
limit often used in the literature is valid for the electron 
reservoirs having a relatively small tunneling rate (negli- 
gible leakage effect) and a large spectral width (negligible 
memory effect). The former implies the validity of the 
Born approximation and the latter corresponds to the 
Markovian approximation. The chemical potentials of 
the electron reservoirs controlcd by the external bias volt- 
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FIG. 9: The same plot as in Fig.|S]by varying d^^R but fixing 
Fl.r = A and other parameters are those used in Fig. [S] 



FIG. 10; The lines to separate the Markovian and the non- 
Mar kovian dynamics in the Fn—t plot, each line corresponds 
to a given value of dL,R- The parameters /j. — E — A — lO^eV 
are used. 



age just modifies the values of the dissipation-fluctuation 
coefficients without altering the characteristic times 
of the reservoirs and the system, as shown in Figs. [3] and 
[71 However, the chemical potential can be very efficient 
to suppress the leakage effect. Thus, /i — _B is a competi- 
tive control parameter in the coherence control of charge 
qubit, as we will see later. 

In order to show clearly when non-Markovian or 
Markovian processes play a major role in charge qubit 
decoherence, we take Fn as an example to examine at 
what time this dissipation-fluctuation coefficient reaches 
its steady value by varying dL,R and F^jj. The result 
is plotted in Fig. [TOl The lines signify transition times 
between the time-dependent fluctuating and the steady 
dissipation-fluctuation coefficients by varying tunneling 
rate at a given spectral width. Non-Markovian dynam- 
ics can be seen mostly in the time range under the lines. 
It shows that the strong non-Markovian dynamics corre- 
sponds to a relatively small spectral widths d^ ji (a strong 
memory effect) and a relatively large tunneling rates Tl h 
(a large leakage effect) compared with the inter-dot tun- 
nel coupling A. Non-Markovian dynamics disappears for 
a large dL,R (> 2A) and a small Tl^r (< A/2). In the 
parameter range of interests to the experiments [ij, the 
non-Markovian processes do not go over more than five 
Rabi cycles, with which a significant effect can be seen 
in maintaining charge coherence, as we will see below. 



B. Decoherence dynamics of charge qubit 

Have examined the time dependencies of all the trans- 
port coefficients (the shifted energy level splitting e'{t), 
the renormalized inter-dot tunnel coupling A'(t) and 
the dissipation-fluctuation coefficients F(i),F'^(i)) in the 
master equation for both biased and unbiased double 



dot, we shall discuss now the decoherent dynamics of 
the charge qubit in this subsection. Experimentally the 
coherence manipulation of charge qubit is performed for 
the double dot on resonance, e = Ei — E2 = 0- The 
corresponding shifted energy level splitting £'{t) remains 
zero. Then the energy eigenbasis of the charge qubit 
refers actually to the molecular anti-bonding and bond- 
ing states, namely, |±) = -^(|1) ± |2)). The oscillation 
between the coherently coupled localized charge states 1 1) 
and 1 2) as coherent superpositions of the molecular states 
describes the charge coherence, where the renormalized 
Rabi frequency ^l{t) — A'(t) is just the shifted inter-dot 
tunnel coupling for the symmetric double dot. The time- 
dependent dissipation-fluctuation coefficients T{t),T'^{t) 
will disturb this coherent oscillation and cause charge 
qubit decoherence. 

To be specific, we let the initial state be pn = 1 and 
examine the time evolution of the density matrix under 
various conditions. We calculate first the rate equations 
([55)1 for the no-Coulomb-interacting double dot. The typ- 
ical population evolution shown in Fig. [TT] tells that the 
double occupancy is favored. In fact, the charge qubit of 
a double dot is designed in the strong inter-dot Coulomb 
blockade regime where the state of simultaneous occupa- 
tion of two dots is excluded. In other words, unlike the 
charge qubit initialization where both the interaction-free 
and the strong interaction rate equations give qualita- 
tively the same result, for the coherence control of the 
charge qubit the interaction-free rate equations (|38p is 
invalid. The charge qubit dynamics must be described 
by the strong-interaction rate equations ([47|) . For the 
rate equations P7)) in the strong inter-dot Coulomb re- 
pulsion regime to be held for charge qubit manipulation, 
the energy difference between the fermi surfaces of the 
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FIG. 11: The time evolutions of the populations pii(the red 
solid line), p22(the blue long-dashed line), poo (the black short 
dashed line) , and P33 (the green dashed-dotted line) calculated 
using the interaction-free rate equations (|38p in the symmet- 
ric double dot with A = lOpeV, fj.- E = 30peV, dL,R = A/2 
and Fl.r = dL,R- Without considering the inter-dot Coulomb 
repulsion, raising up of the chemical potentials ^ > E accu- 
mulates charges into the double dot. 



reservoirs and the energy levels of the dots, fi — E, can- 
not be too large. The inter-dot Coulomb repulsion in 
the samples is estimated to be 200/ueV The value 
of /i — £' that can be taken most safely should be not 
larger than 50/xeV, a quarter of the inter-dot Coulomb 
repulsion energy. If — i? is taken over 100 neV, it is 
comparable to the Coulomb repulsion energy so that the 
doubly occupied state cannot be completely excluded. 
For convenience and consistency with the discussion in 
the previous subsection, we still take the quantum dot 
parameters Ei = E2 = E and A = lO/xeV. As one will 
see with fi — E = 25 ^ 50/^eV, and Tl,r and dL,R be- 
ing in a reasonable range, the charge qubit can maintain 
coherence very well. 

In Fig. [T^l we plot the time evolution of the reduced 
density matrix by varying the aligned fermi surfaces. The 
result shows that the coherent oscillation of the charge 
qubit depends sensitively on the hight of aligned fermi 
surfaces from the resonant levels of the double dot, i.e., 
ji — E. When ji — Eis not too large (< A) and r^^ij is not 
too small (> A), the population pn decays very quickly. 
This is because although the state of simultaneous oc- 
cupation of two dots is excluded, there is still chance 
for electrons to escape from the dots into the reservoirs 
such that both dots become empty, namely poo 7^ [sec 
Fig. fTWb)] as a leakage effect. This effect can be sup- 
pressed when the fermi surfaces arc aligned such that 
fi ~ E must be relatively larger than r^ jj. As a result, 
the charge qubit can maintain the coherence very well. 
In Fig. [TWa). we see that the charge coherence is per- 
fectly maintained for ~ E = 50^eV. Meanwhile, the 
real part and the imaginary part of the off-diagonal den- 
sity matrix element pi2 exhibit quite different dynamics. 
The imaginary part of pi2 has a similar oscillatory fea- 
ture as pii [see Fig. [T^ d)]. which depicts the coherent 



tunnel coupling between two dots. While the real part 
of P12 goes down to be negative [see Fig. [r^c)] which is 
related to the loss of the energy (dissipation). We find 
that to maintain a good coherent dynamics for charge 
qubit, the fermi surfaces of the reservoirs is better to be 
aligned above the energy levels of the double dot not less 
than 2 A. 
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FIG. 12: The time evolutions of the density matrix at different 
fi ~ E with E fixed. The red solid line is for fi — E = 0, the 
blue long-dashed line is for fj. — E = lOpeV and the black 
short-dashed line is for jj, ~ E — 50peV. dL,R = A/2 and 
Tl,r = A. 



In Fig. [13] we plot the time evolution of the reduced 
density matrix at a few different tunneling rates but 
fixing other parameters. Experimentally, the tunneling 
rates between the dots and the reservoirs are also tunable. 
We fixed the spectral widthes at cIl.r = 25//eV for which 
the memory effect is largely suppressed. Meantime, we 
take ji — E = 50/icV (much larger than Vl^r) such that 
the leakage effect is also largely suppressed. When the 
tunneling rates are small (< A/2), the time-dependence 
of these transport coefficients in the master equation are 
negligible (as shown in Fig. [5]) so that no non-Markovian 
dynamics can be observed. The corresponding charge 
dynamics is given by the red solid lines in Fig. [T3] where 
the oscillation frequency is time independence, consis- 
tent with the result in Markovian approximation. When 
the tunneling rates become large (> A), the charge fre- 
quency is largely shifted and varies in time. The larger 
the tunneling rates are, the faster the electron oscillates 
between two dots and the reservoirs, thus the stronger 
the non-Markovian dynamics occurs. Here the decay of 
the coherent charge oscillation (see the blue long-dashed 
and black short-dashed lines in Fig. [T5|) is not due to 
the charge leakage (which has been mainly suppressed 
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by raising up the fermi surfaces) but a back-action deco- 
herence effect of the reservoirs when the tunnehng rate 
becomes larger. Increasing the tunnehng rates leads to 
more charge leakage. However, comparing the magni- 
tude of Poo with that of pn shows that the damping of 
coherent charge oscillation in the presence of a large tun- 
neling rate is not primarily due to charge leakage but a 
non-Mar kovian back-action decoherencc effect. 
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FIG. 13: The time evolutions of the density matrix at different 
coupling constants. The red solid line is for Tl,r ~ 2.5/ieV, 
the blue long-dashed line is for Tl,r = WfieV and the black 
short-dashed line is for Tl,r ~ 25/xeV. The chemical poten- 
tials are kept a,t fi — E — 50/ieV and (Il.r = 25/ieV. 



The more non-Markovian dynamics can be seen by 
varying the spectral widths cIl.r- When the spectral 
width is comparable to the inter-dot coupling A, namely 
the characteristic time of reservoirs is comparable to the 
characteristic time of the double dot, the non-Markovian 
dynamics becomes the most significant in the time evo- 
lution of charge coherence. Although it is currently not 
clear how to tune the spectral widths in experiments, it 
is still interesting to see what roles the spectral widths 
(or more generally speaking, a non-constant spectral den- 
sity) play in the charge decoherence dynamics. We plot 
in Fig. [13] the time evolutions of the density matrix el- 
ements at various spectral widths with fi — E = bOficV 
and Fl,/? = A. As we can see if d^ji is small (< A/2), 
the coherent charge dynamics is well preserved although 
it is a strong non-Markovian process. With increasing 
the spectral widths, the decoherent charge dynamics be- 
comes visible and also becomes Markov type. Widening 
the spectral widthes damps the coherent charge oscilla- 
tion. Wider spectral width also causes more charge leak- 
age. But comparing the magnitudes of poo with pu in 



Fig. [M] shows again that the damping of coherent charge 
oscillation in the presence of a wide spectral width is a 
Mar kovian back-action decoherence effect. The frequen- 
cies of coherent charge oscillation are also shifted differ- 
ently from the bare Rabi frequency for different spectral 
widths presented, [see Fig. [UTal-fd')]. 




FIG. 14: The time evolutions of the density matrix at different 
spectral widths. The red solid line is for dL,R ~ IpeV, the 
blue long-dashed line is for dL,R ~ 5fieV and the black short- 
dashed line is for dL,R ~ 25peV. p — E is kept at 50peV and 
Tl.a = A. 



From the above analysis, we find that the charge qubit 
coherence can be maintained very well when either the 
spectral widths dL,R or the tunneling rates Tl.r are suf- 
ficiently smaller than p — E. However, when both the 
spectral widths dL^R and the tunneling rates Tl,r. be- 
come comparable to p — E, the decay of charge coherency 
can be seen within a few cycles of the bare Rabi oscil- 
lation. The difference between the non-Markovian and 
the Markovian process manifests in the difference be- 
tween the renormalized Rabi frequency U,{t) = A'(t) (for 
the symmetric double dot) and the bare Rabi frequency 
fio = In Fig- [151 we plot the average renormalized 
Rabi freuency (57) = Viit) by varying the spectral widths 
dL,R and the tunneling rate T l r, respectively, and make 
a comparison with the bare Rabi frequency (the black 
short dashed lines). It shows that the renormalized Rabi 
frequency has a large shift from the bare one, except 
for the region (the Markovian regime) with the spectral 
width dL,R > 2A and the tunneling rate F^.^^ < A/2, 
where the renormalized Rabi frequency is close to the 
bare one. In Fig. [TH we plot the time evolution of den- 
sity matrix pn using the strong-interaction rate equa- 
tions ([TZl) for a few different sets of L^R^d^^R) and 
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make a comparison with the Markovian approximation 
determined by the rate equations (gS]) , where the charge 
coherence dynamics in the double dot from Markovian to 
the non-Markovian processes is clearly demonstrated. 
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FIG. 15: The average renormalized Rabi frequency (Q) — 
Q{t) by varying the spectral widths cIl.r and the tunneling 
rate Tl,r, respectively, and comparing with the bare Rabi 
frequency (the black short dashed lines). Here we take fi—E = 
50/ieV to suppress the charge leakage. 
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FIG. 16: The time evolutions of pn obtained from the rate 
equation (|47|) and the corresponding Markovian limit (|48|) for 
the strong-Coulomb-interacting double dot. (a) the strongest 
non-Markovian regime for a small spectral width and a large 
tunneling rate, (b) the non-Markovian dynamics controlled 
by a large tunneling rate alone, (c) the non-Markovian dy- 
namics controlled by a small spectral width alone, (d) The 
Markovian regime corresponding to a small tunneling rate and 
large spectral width. Here we take jj. — E = SOfieV. 



C. Relaxation time Ti and dephasing time T2 

To understand quantitatively the decoherence dynam- 
ics of charge qubit in the double dot, we shall now ex- 
tract the decay rates (the relaxation time Ti and the 
dephasing time T2) for various manipulation conditions. 
The relaxation time Ti characterizes the time going from 
the anti-bonding state (with higher energy) to the bond- 
ing state (with lower energy) in the energy eigenbasis 
|±) = ^(|1) ± |2)) (for symmetric double dot). It 

is descritaed by the decay of the diagonal element of 
the density matrix in the energy eigenbasis, = 
■^(Pii + P22) + R.c(pi2). In the ideal case (such as in 
NMR), pii +P22 = 1 so that Ti is completely determined 
by the real part of pi2. For the charge qubit in double 
dots, in general pn + P22 < 1 because of charge leakage. 
But in the practical manipulation, the fermi surfaces of 
the reservoirs are set high enough from the resonant en- 
ergy levels of the double dot such that charge leakage can 
be suppressed. Thus we can still extract Ti from Repi2. 
The decoherence (dephasing) time T2 corresponds to the 
decay of the off-diagonal element of the reduced density 
matrix = ^{Pu - P22) - *Im(pi2)- As we will 

see later the time-dependencies of pii,p22 and the imag- 
inary part of P12 behave very similar. This may tell us 
that T2 can be extracted either from pii,p22 or Impi2. 
Experimentally, one extracted T2 from P22 by measuring 
the current proportion to P22 01 ■ 

Having made the above analysis, we shall extract the 
decay rates of the charge coherent oscillations from the 
time evolution of the reduced density matrix elements 
Kepi2{t) and lmpi2{t) or pii(t) by fitting a decay oscil- 
lating function plus an offset to the numerical data. An 
intuitive fitting function for the charge oscillation decay 
would be Ae~^* cos((f2)i) -|-C. However this fitting func- 
tion fails for pii{t). The typical damping oscillation of 
Pii{t) shows that it converges to a steady value of pn < i 
at large t. Thus pu {t) can be well described by the fitting 
function 

PiiW = hit) cos^ i{n)t/2) + flit) sm\{n)t/2) (52) 

where ,fh,iit) = Ah,i expi-Bh,it^) + Ch,i are used to fit 
the downward shift of the peaks and the upward shift 
of the valleys in the damped oscillation, respectively. 
The oscillating function is the squares of sine and cosine 
functions with half the oscillation frequency, sin^((ri)t/2) 
and cos^((r2)</2), rather than cos((ri)i) and sin((ri)t). 
This can be easily understood by considering an ideal 
qubit. Its Rabi oscillation conditioned to the initial state 
. f,s fcos^nt/2) -|sin(m)\ 
pn(0) = 1 IS given by pit) = , ^.^^^^^^ sil\nt/2) ) 

in the localized charge state basis. Similarly, the off- 
diagonal density matrix elements can be described by 
the fitting functions 

Impi2(t) = {A, exp(-Sjt") + C^} sin((f7)0, (53a) 
Repuit) = Ar expi-Brt') + Cr. (53b) 
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The fitting parameters A^, Bx,Cx for x = h,l,i,r are 
generally different for different fitting function and dif- 
ferent data. Fig. [IT] is a plot of using the fitting func- 
tion of ([5^ to fit the exact numerical solution of pii{t). 
The results show that with the decay function fh{t) = 
0.48 exp (-0.58t) + 0.50 and fi{t) = -0.47exp(-0.54t) + 
0.49, the fitting function gives very much the same 
solution as that obtained numerically from Eq. (|47p for 




FIG. 17: (a) Using tlie fitting function fh{t) cos^{{n)t/2) + 
fi{t) sin^((f2)t/2) (the red solid line) to fit the exact numerical 
solution of p\i{t) (the green dashed-dotted line). Here we use 
the parameters ^ — E — 50/ieV, Tl,r = 25/xeV and cIl.r ~ 
25MeV. 



In order to have a better fitting to various tunable pa- 
rameters, we plot in Fig. [TH] the fitting errors for pii(t) 
by averaging the deviations between fh (t) and the peaks 
of the oscillating pii{t) in a large range of the chemi- 
cal potential, the tunneling rate as well as the spectral 
width. The top two plots in Figs. HH] are the fitting er- 
rors by varying pl^r but fixing (Tl^h, dL,R) = (25, 5)//eV 
(the strong non-Markovian regime) and (5, 25) /icV (the 
Markovian limit regime). We find that in the strong 
non-Markovian regime, when the fermi surfaces are not 
far away from the resonant levels of the double dot 
[ji — E < 25/ieV), the best fitting is a sub-exponential 
decay with s < 1. When p — E > 25/ieV, the fitting 
function becomes a simple exponential decay function 
{s — 1), while, in the Markovian limit, the best fitting is 
just a simple exponential decay for all the values of fi—E. 
The middle two plots in Figs. [TH] are the fitting errors by 
varying Tl,r but fixing {p - E, cLl^r) = (10, 5)^eV (the 
non-Markovian regime) and (30, 25) peV (the Marko- 
vian limit). Again we see that in the non-Markovian 
regime, the best fitting is a sub-exponential decay with 
s < 1 for all the values of T^ ji except for some very small 
^L,R (< ^/2) which indeed enters the Markovian limit 
where the fitting function becomes a simple exponential 
decay function. In the Markovian limit (the large spec- 
tral width limit here) , the best fitting is given by a simple 
exponential decay (s = 1) again. The bottom two plots 
in Figs. [15] are the fitting errors by varying d^ ji but fixed 
{p — E,T L^ft) = (10, 25)/ieV (the non-Markovian regime) 
and (30, 5) /icV (the Markovian limit). It tells that in 



the non-Markovian regime, when cLl^r < A (the strong 
non-Markovian regime), the best fitting is still given by 
a sub-exponential decay (s < 1), while for d^^R. > A the 
system transits to the Markovian, the fitting function be- 
comes again a simple exponential decay function. In the 
Markovian limit (small tunneling rate limit), the best fit- 
ting is just given by a simple exponential decay (s = 1), 
as one expected. 
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FIG. 18: The average fitting errors by fitting the exact nu- 
merical solutions of pii(t) with the fitting function given by 
Eq. (|52p . The top two plots are the fitting errors by vary- 
ing jj. — E but fixing {T L,R,dL,R) = (25, 5)/ieV and (5, 25) 
peM , respectively; The middle two plots for varying Tl,r 
but fixing — E,dL,R) = (10, 5)/ieV and (30, 25) /leV; and 
the bottom two plots are obtained by varying dL,R but fixed 
{^-E,Vl,r) = (10, 25)/ieV and (30, 5) peV . The parameters 
are chosen such that the left three figures correspond to the 
non-Markovian regime while the right three figures are the 
Markovian limit, s = 1/2 corresponds to the sub-exponential 
fitting (red lines), s = 1 is the simple exponential fitting (blue 
long-dashed lines), and s = 3/2 is the super-exponential fit- 
ting (black dashed lines). 



The above analysis shows that Markovian decoher- 
ence processes lead to an exponential decay and a sub- 
exponential decay seems to occur mainly in strong non- 
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Markovian processes. But this does not imply a sim- 
ple exponential decay being necessarily Markovian. In 
Fig. (fT9|) , we plot the average fitting errors of the fitting 
functions (j53p with the exact numerical solution of the 
real and imaginary parts of the off-diagonal density ma- 
trix element pi2(t) in the strong non-Markovian regime 
with the spectral widths cIl^r ~ 5/ieV and the tunnel- 
ing rates Tl ji = 25^eV. The results show that for both 
Impi2(<) and Repi2(i), the best fitting is given by simple 
exponential decay for the whole range of chemical poten- 
tial up to fj. — E = 40^eV. This forces us to carefully 
look at the results in Fig. [TH] We find that all the results 
with a sub-exponential decay show up for small ^ — E 
values (< 2 A) where the charge leakage effect cannot be 
effectively suppressed by the chemical potentials of the 
reservoirs. This tells that the sub-exponential decay in 
pii(i) is a charge leakage effect. The slightly different de- 
cay behaviors for pii{t) and lmpi2{t) actually come from 
the charge leakage effect contained in pii(i) rather than 
a consequence of non-Markovian dynamics, as we will see 
more later. 
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FIG. 19: The average fitting errors by fitting the exact nu- 
merical solution of Impi2(t) and Repi2(t) with the fitting 
functions given by Eq. (|53p . where we vary — E but fixing 
dL,R = 5^eV and Vl,r ~ 25/ieV. The red lines correspond 
to the sub-exponential fitting function (s = 1/2), the blue 
long-dashed lines are the simple exponential fitting function 
(s — 1), and the black dashed lines are the super-exponential 
fitting function (s = 3/2). 



Now we shall extract the relaxation time Ti from 
Re/5i2(t) and the decoherence time T2 from Impi2(t) or 
Pii{t) using the concept of half life from the exact nu- 
merical solutions and from the fitting functions (|53[) and 
(PI). The results is plotted in the Figs. [2D]1211 In Fig.HD] 
we plot the relaxation time Ti and the decoherence time 
T2 from Repi2(t) and Impi2(t); respectively, by varying 
the chemical potential fi — E, the tunneling rates F^ jj 
and the spectral widths (Il^r differently. Fig. [^nra)-fb) 
are the plots of Ti and T2 by varying p — E but fixing 
i^L.RjdL.R) ~ {25,5)peV (solid fines, corresponding to 
the non-Markovian processes) and (5, 25) peV (dashed 
lines for Markovian limit). The results tell that the deco- 
herence effect (for both Ti and T2) is large for the small 
chemical potential p — E due to the large charge leak- 



age effect. Increasing p — E reduces the charge leakage 
effect, thus also reduces the decoherence effect. When 
p — E IS larger than Vl,r and d^ R [goes up to 30 pcN 
in Fig. [2nKa)-(b)], the decoherence effect quickly reaches 
to a minimum value (the longest decoherence time ^ 2 
ns). Fig. Wic)-{d) plot Ti and T2 by varying Tl,r but 
fixing {p — E, (Il^r) = (10, 5) (solid lines, where both the 
charge leakage effect and the memory effect are supposed 
to play an important role) and (30, 25) peV (dashed lines, 
where both the charge leakage effect and the memory ef- 
fect arc ignorable). It shows that for a small tunneling 
rate between the reservoir and dot (< A) the decoher- 
ence effect is weak. Increasing T^^r enhances the non- 
Markovian dynamics effect and also enhances the deco- 
herence (shorting the relaxation and dephasing times). 
Fig. EOf e)-ff) plot Ti and T2 by varying d^ R but fixing 
{p — E, Tl r) — (10, 25) (solid lines, where the charge 
leakage effect is dominated) and (30, 5) pcV (dashed 
lines, where the charge leakage effect is negligible). We 
find that when the memory effect must be considered 
(corresponding to a small spectral width), the decoher- 
ence effect is small (or the relaxation and dephasing times 
become longer). Increasing dL^R reduces the memory or 
non-Markovian dynamics effect and enhances the deco- 
herence (shorting the relaxation and dephasing times). 
It is interesting to see that in all the cases, the relaxation 
time is close to the dephsing time: Ti ~ T2 in the Marko- 
vian limit (dashed lines), where for the non-Markovian 
regime, the relaxation time is larger than the dephasing 
time up to twice: Ti < 2T2 (solid lines). 

As we have pointed out that the decay behaviors for pn 
and Im/3i2 are slightly different from the fitting function 
in the non-Markovian regime. In Fig. 1211 we compare the 
decoherence time T2 extracted from pii(t) and Impi2(t), 
respectively, with the same conditions as used in Fig. 1201 
namely by varying the chemical potential p — E, the tun- 
nefing rates Tl,r and the spectral widths dL.R differently. 
As we see although quantitatively the dephasing time ex- 
tracted from pii(t) and Impi2(i) are in the same order, 
but there are some obvious differences in certain range of 
the tunable parameters where the charge leakage effect 
play an important role. This is indeed clearly show in the 
right-bottom two plots in Fig. [21] where p — E = 30peV. 
When Tl^r are small (< A/2) so that the charge leakage 
is negligible, the dephasing time T2 extracted from pii(t) 
and Impi2 (t) are very close to each other over there. This 
again indicates that it is the charge leakage effect that re- 
sults in the slightly different decay law for pii{t) (follows 
a sub-exponential decay) and lmpi2{t) (by a simple ex- 
ponential decay) in the non-Markovian regime. 

In fact, the original definition of the decoherence for a 
qubit is given by the decay of |(+|p|— )| though it is not 
easy to be measured directly in experiments. Ultimately 
when the charge qubit is completely decohered, pu ~ P22 
and Im(pi2) ~ thus |(+|p|— )| at the asymptotic 
time. We have verified this property in our exact numer- 
ical calculation. Thus the off-diagonal reduced density 
matrix clement in the energy eigenbasis. |(+|p(OI~)L *3an 
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FIG. 20: The relcixation time Ti and the dephasing time T2 
extracted from Repi2(t) and Impi2(t), respectively, (a)-(b) by 
varying the chemical potential fx — E but fixing {fi — E, (Il^r) ~ 
(10,5) (solid lines) and (30, 25) ^eV (dashed lines), (c)-(d) 
by varying the tunneling rates Tl,r at {fj, — E,dL,R) = (10, 5) 
(solid lines) and (30, 25) fj,eV (dashed lines), and (e)-(f) by 
varying the spectral widths dL,R but fixing (/i — E,Tl,r) = 
(10,25) (solid lines) and (30, 5) /xeV (dashed lines). 



FIG. 21: The decoherence time T2 extracted from pii(t) 
and lmpi2{t), respectively, with the same conditions used in 
Fig. 1201 namely by varying the chemical potential /j. — E, the 
tunneling rates Vl.r and the spectral widths dL,R differently. 



be well fitted by Ae-^*° + C with C = 0. In Fig. we 
compare the results of T2 extracted from )| and 

lm{pi2){t), respectively. It is remarkable that the de- 
phasing time T2 obtained from )| and lin{pi2){t) 
are almost exactly the same in a wide range of parame- 
ters concerned here( fj, — E is from to 40/xeV, Tl^r and 
dL,R are from 1 to 30/ueV at A = 10/xeV). The decoher- 
ence (dephasing) time T2 obtained here is between 0.2 to 
2 ns, except for the case shown in the left-bottom plot in 
Fig. 1211 where the decoherence time even smaller. 

Now we shall end this section with a brief summary 
in the following. The non-Markovian coherence and de- 
coherence dynamics of charge qubit is dominated by two 
major effects, the memory effect and the leakage effect 
in the double dot gated by electrode reservoirs. The for- 
mer becomes a dominate effect when the time scale of 
the reservoirs is comparable to the time scale of the dou- 
ble dot. The latter becomes an important effect when 
the electron tunneling strength between the reservoirs 
and dots is tuned to be large. These two characters are 
suitably described by the spectral widths and the tun- 



neling rates embedded in the Lorcntzian spectral density 
([25|l we used. Strengthening the couplings between the 
reservoirs and dots and widening the spectral widths of 
the reservoirs disturb the charge coherence in the dou- 
ble dot significantly. However, reasonably raising up the 
chemical potentials pl.r can suppress charge leakage and 
maintain charge coherence. The left uncontrollable de- 
coherence factor is the spectral width which characters 
how many electron states in the reservoirs effectively in- 
volving in the tunneling processes between the dots and 
the reservoirs. The smaller the spectral width is (the less 
the electron states involve in the electron tunneling) , the 
better the charge coherence can be maintained. The de- 
cay of the charge coherent oscillation is well described by 
a simple exponential decay for the off-diagonal reduced 
density matrix elements, but the diagonal ones (popu- 
lations) are better described by a sub-exponential law 
when charge leakage is not negligible. Otherwise simple 
exponential decay is better for both the non-Markovian 
and Markovian regimes. The relaxation time Ti and the 
dephasing time T2 can be extracted from the exact nu- 
merical solution Rcpi2(t) and Impi2(t); respectively, with 
the result Ti < 2T2 for a broad parameter range we used. 
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FIG. 22: The decoherence time T2 extracted from |p+,_(f)| — 
|(+|p|— )| and Impi2(t), respectively, with the same conditions 
used in Fig. 1211 namely by varying the chemical potential 
fi — E, the tunneling rates Tl.r and the spectral widths cIl.r 
differently. 



V. CONCLUSION AND DISCUSSION 

In this paper, we have developed a non-perturbation 
theory to describe decoherence dynamics of electron 
charges in the double quantum dot gated by electrodes. 
We extended the Feynman- Vernon influence functional 
theory to fcrmionic environments and derived an exact 
master equation ([33l) for the reduced density matrix of 
the double dot without including the inter-dot Coulomb 
repulsion at beginning. The contributions of quantum 
and thermal fluctuations induced by the electron reser- 
voirs are embedded into the time-dependent transport 
coefficients and ([57]) in the master equation. These 
time-dependent transport coefficients are completely de- 
termined by the non-perturbation dissipation-fluctuation 
equations of motion (|30[) . The exact master equation is 
then further extended to the double dot in the strong 
inter-dot Coulomb interacting regime in terms of Bloch- 
type rate equations (|17|) where the strong Coulomb repul- 
sion simply leads one to exclude the states correspond- 
ing to a simultaneous occupation of the two dots from 
Eq. ([55)) . Our theory is developed for a general spec- 
tral density of the reservoirs at arbitrary temperatures 



and bias. Other approximated master equations used 
for the double quantum dot can be obtained at well de- 
fined limits of the present theory. This non-perturbation 
decoherence theory allows us to exploit the quantum de- 
coherence dynamics of the charge qubit brought up by 
the tunneling processes between the reservoirs and dots 
through qubit manipulations. 

We then used the master equation (in terms of the 
rate equations) to study the non-Markovian decoherence 
dynamics of the double dot charge qubit with the back- 
action of the reservoirs being fully taken into account. To 
make qualitative and also quantitative understandings 
of the charge qubit decoherence, we numerically solve 
the dissipation-fluctuation integro-diffcrcntial equations 
of motion using a Lorentzian spectral density. We exam- 
ine the time dependence of all the transport coefficients 
from which the time scales within which non-Markovian 
processes become important in the charge coherent dy- 
namics are determined. The correlation time of the elec- 
tron reservoirs (in terms of the spectral widths in the 
Lorentzian spectral densities) and the electron tunnel- 
ing strengths between the reservoirs and dots (in terms 
of the tunneling rates in the Lorentzian spectral den- 
sities) characterize the time scales for the occurrence 
of non-Markovian processes. Non-Markovian processes 
dominate the charge coherent dynamics when the spec- 
tral width is comparable to the inter-dot tunnel coupling 
where the memory effect plays an important role, and/or 
when the tunneling rates between the reservoirs and dots 
become strong such that charge leakage becomes a main 
effect for decoherence. Raising up the fermi surfaces of 
the reservoirs can suppress charge leakage. The Marko- 
vian limit can be reached with a weak tunneling rate 
and a large spectral width, where perturbation theory 
becomes valid and the spectral density is reduced to a 
constant. The decay of the charge coherent oscillation 
is well described by a simple exponential law, except for 
some special regime where the charge leakage is not negli- 
gible such that the evolution of state populations is better 
described by a sub-exponential decay. We also extracted 
the relaxation time Ti and the decoherence time T2 con- 
sistently from different elements of the reduced density 
matrix and obtained a general result Ti < 2T2 which is 
of the order of ns or less in a broad parameter range we 
considered. These results arc ready to be examined in 
experiments 

Although we concentrate in this paper the electron 
charge coherence (decoherence) dynamics in the double 
quantum dot system, the theory we developed in this 
work can also be applied to investigate other physical 
properties, such as quantum transport phenomena, in 
various quantum dot structures. In fact, the spectral 
density describing the spectral distribution of electron 
reservoirs and electron tunneling processes between the 
reservoirs and dots has not been well determined exper- 
imentally. We used a Lorentzian spectral density that 
has been used by others yet still been waiting for justi- 
fying in experiments. The tunneling rates T^n in the 
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Lorentzian spectral density is indeed tunable in experi- 
ment, and have been extracted from current spectra by 
assuming a constant spectral density. With the capabil- 
ity of monitoring the time evolution of electronic popula- 
tion transfer, we can look at closely the short time trans- 
port properties in this double dot device from which we 
may extract the tunneling rates and the spectral widths 
for the Lorentzian spectral density or other possible forms 
of the spectral density. Otherwise it may be other sur- 
roundings (phonons and fluctuation in impurity configu- 
rations, etc.) that play an important role in the dynamics 
of charge qubit decoherence. These together with other 
transport properties in various nanostructures, such as 
Kondo effect and Fano resonance, etc., deserve a sepa- 
rate study. We will leave these research problems to be 
addressed in separate papers. 
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APPENDIX A: DERIVATION OF THE 
INFLUENCE FUNCTIONAL 



Let the electron reservoirs be initially in a thermal 
equilibrium state, then 



(/o|Pi=;(^o)l3o> =n<^"l7e"''^'""''^"''=""l.9'^o) (A3) 



Ik 



where fJ-L.R are the chemical potentials of the source and 
drain electron reservoirs connected to the dots 1 and 2, 
respectively, and Z is the fermion partition function of 
the reservoirs Z = Hifcle"^^'^"""^'^ + !)■ Obviously, the 
influence functional can be written as J^[£.*,S,,'>J*,'ri\ = 
riife -^ife [C* ) ?j ^7* I ^] • Furthermore, since the Hamiltoni- 
ans are quadratic, the path integrals in the influence func- 
tional can be exactly calculated using either the Gaussian 
integrals or the stationary path method. Here we present 
the clculation based the stationary path method. The 
forward stationary pathes of the electrons in the reser- 
voirs is determined by 



The propagating function governing the time evolution 
of the reduced density matrix is given by Eq. (|22p . in 
which the generalized Feynman- Vernon's influence func- 
tional is defined by 



+ S'7[r,«,/*J]-S;[r,*,7),g*,s]}^ (-^^^ 

where pE is the initial density operator of the fermion 
reservoirs, ./o,5o,/Af and their complex conjugates are 
Grassmann numbers introduced in the fermion coher- 
ent state representation, Se is the action of the electron 
reservoirs, and Sj stands for the action of the interaction 
between the dots and the reservoirs. Explicitly, 

f;kit)fik{t) + f;kito)fik{to) 



^SE[rJ]= { 

l=l,2;k 

' flkflk ~ flkflk 



dr 



to 



i^lk flkflk 



] 



~ ^y^.sE,ik[fiki fik], 

Ik 

Slid, f\f] = -Y.f driUik^Uik + tlJtk^i) 
Ik "'to 



(A2a) 



SlMG^^i^ flk^ flk] 



(A2b) 



fikir) + ieikfikir) = -it-zfcCrW, 
fik{T-) - ieikftkij) = iUikCi(j), 



(A4a) 
(A4b) 



with i = 1, 2 for / = L, i?, respectively. The solutions to 
the stationary path equations (|A4p are 



fikir) = fik{to)e-'"'^^-">^ - / dr'e— (---')C.(r') 



f^kir) = fUty'"-^^-'^ - ^tak / dr'e--(^"^')C(r') 



(A5) 



With the similar solutions for the backward stationary 
pathes, the Ik component of the influence functional is 
then given by 



Ik 



24 



dT 



dr' (e-- (---')77* (r')??, (r) + ^^-^'^^ (r)^, (r')) 



+ / rfr'e---(---')(,7*(rK.(r')-/(£zfc)K(T)+77*(r)]K.(r')+r7,(r')])]}, (A6) 



r 



where f{sik) is the fcrmi distribution function, f{eik) = 
Of,,, "'",,,1 , , . Sum up contributions from aU fermion 

e'^^ Ik i^u 

modes in the electron reservoirs, we get the influence 
functional Eq. (j23p with the dissipation-fluctuation ker- 
nels given by Eq. 



APPENDIX B: DERIVATION EOF THE EXACT 
PROPAGATING FUNCTION 

In this section we show how to use the solutions of the 
equations of motion, (|28p . to determine the time depen- 



dent coefficients in the master equation. The equation 
of motion for 2 is just the complex conjugate equa- 
tion of 7712 while that of rjl 2 is the complex conjugate to 
the equation of ^1^2- Meantime both ^*2(''') ^^nd ?7i.2(''") 
are fixed at r = t and both j]l 2 ("J") a-nd ^i,2(''') are fixed 
at r = to- We only need to solve the set of equations of 
motion for ^1.2 (t) and 771.2 (''')• The solutions to the equa- 
tions of motion for r]l 2{''') and Ci2(''') "-an be obtained 
by conjugating the solutions of ^i,2{t) and 771,2 (t) with 
corresponding replacement of the boundary conditions. 

Let x(''") =£,{''')+ vi'T')^ the equations of motion 
becomes 



J 



X(r) + M t I; ) X(r) 



to 



, / (Fil{t~t') \ 







F2r{ 



(Bla) 
(Bib) 



r 



where we have used the following notations for brevity, Eq. (|B1[) can then be expressed in terms of the new vari- 
ables m(t), u{t) and v{t) as 



6(r)j' U(r) 



(B2) 



To solve the above equations of motion, we introduce the 
new variables u{t)^ ''^{t) and v{t) such that 

x{r)^u{T)x{t), (B3a) 
i{T)=u{T)i{to)+v{T)x{t). (B3b) 







to 



(Fil{t-t') u , , 

' F2R{T-T')y^^> 



to 



F^dr-r') 





(t- 



(B4a) 
(B4b) 

(B4c) 
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with the boundary conditions Uij{t) = 5ij, Uij(to) = Sij 
and Vij{to) = 0, respectively. Obviously, Eq. (|B4a[) is 
the backward version of Eq. ([Bibl) . Therefore, u{t) = 
u^{t + to - t) for to<T <t. 

Now Ci,2(i), '?r,2(*o), ?7i,2(to) and ?7^_2(0 can be fac- 
torized from the the boundary conditions, £,{to) = ^o, 
ri{t) = j]f. Explicitly, let t = to for Eq. (|B3aP and t = t 
for (|B3bp . we have 

Co + v(to) ^ u{hmt) + Vfl (B5a) 
m^umo + v{t)m+rif]- (B5b) 

Note that u{to) = u^{t), the above algebraic equa- 



tion gives the solution ((29|) . Similarly, let ^*(t) = 
(a(^) ^2(^)) and rj*ir) = (771* (r) if,{T))', we have 

- 77* [/ - ut(t)(/ - vHt))-'u{t)] , (B6a) 
7f (t) = [77Syt(i) + ^*.„t(^)] (J _ „t(i))-i. (B6b) 

Substituting the relations (pQ]) and (|B6P into the propa- 
gating function (|27p and using the fact that v(t) is hermi- 
tian at r = i, we obtain the exact propagating function 
((3T|l for the double dot gated by bias electrodes. 
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